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ABSTRACT 



0?his investigation is concerned with the interaction 
between a turbulent flow and certain types of structure re- 
sponding to its excitation. The turbulence is typical of 
those associated with a boundary layer, having a cross-spectral 
density indicative of convection and statistical decay, A 
number of s-'cructural models are considered in the investigation 
Among the one-dimensional models are an unsupported infinite 
beam and a periodically supported infinite beam. The first 
model is used to deVelope the basic ideas which are then 
applied to the more realistic second model resembling 
the fuselage construction of an aircraft. For the two-dimen- 
sional case a simple membrane is used to illustrate the type 
of formulation applicable to most two-dimensional structures. 
However , a small random .variation in the membrane tension is 
included in the analysis since ideally uniform tension never 
exists in practice. Moreover, the mathematical approach used 
in dealing with random membrane tension can be adapted to 
treat other random structural properties in general. Both the 
one-dimensional and two-dimensional structures mentioned above 
are backed by a - cavity filled with an initially quiescent 
fluid to simulate the acoustic environment when the structure 
forms one side of a cabin of a sea- or air-craft. 

It is shown that a decaying turbulence can be con- 


iii 

strucijed ftom superposing infinitely many components, each of 
which is Gonvected as a frozen-pattern at a different velocity, 
Shis superposition scheme reduces greatly the computation time 
by reducing to one-half the number of integration which must 
be performed on a computer, Furthermore , the scheme provides 
a convenient way in which experimentally measured cross- 
spectral density. of the turbulence pressure fluctuation can 
be incorporated directly in the computation, 

The results of the structure -turbulence interaction 
are presented in terms of the spectral densities of the struc- 
tural response and the perturbation Reynolds stress in the 
fluid at the vicinity of the interface. It is found that 
important spectral peaks of the structural response will not 
appear if decays in the turbulence ‘ is neglected in the analysis. 
Thus, the usual Taylor's hypothesis of frozen-pattern turbulence 
is unconservative as far as the assessment of structural 
reliability is concerned. The parirurbation Reynolds stress 
is indicative of the change in the skin- friction drag due to- 
structural motion. It is. shown that, given the statistical 
information of the boundary-layer turbulent pressure field, the 
perturbation Reynolds stress can be changed’ by varying ..the 
structural parameters. Therefore, the present stddy is poten- 
tially useful for designing flight or marine structures to 
minimize the total skin-friction drag. 
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I. IIKPRODUCTIOK 

In recen-fc years, *fcliere has "been considerahle interest 
in the response of panel systems to random pressure fields. 

An aircraft fuselage excited by boundary-layer turbulence or 
by the efflux of a jet engine is a good example. The objec- . 
tives of such investigations Ci-9] vary from minimizing struc- 
tural fatigue damage to reducing noise radiation from the 
panel system to either outside or inside of the cabin. 

Recently, another important application has been suggested, 
namely, to design a panel sjrstem for the exterior of a vehicle 
such that the total s3cin-fr lotion drag force over the vehicle 
is a minimum* 

The typical panel system of an aircraft fuselage is a 
multi-span structure which is characterized by close clustering 
of natural frequencies in well-defined frequency bands. The 
usual normal mode formulation does not lead to practical 
results in this case, since it is almost impossible to calcu- 
late the normal modes of structures with a large number of 
spans due to close proximity of natural frequencies in each 
frequency band. However, if all the panels in a system are 
identically constructed, then the structural configuration is 
spatially periodic and the analysis can be greatly simplified. 
Two alternative methods are available to solve the response 
prcblem of suoh a spatially periodic structures the transfer 
matrix technique S3id the wave-propagation approach 

Cl 7-21]. The computational simplicity of both methods is 
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achieved hy utilizing the fact that the entire system is 
composed of identical subunits. Although these two methods 
are related £ 223 » the transfer matrix techniq.ue is more 
suitable for analyzing a periodic structure with finite total 
length whereas the wave-propagation approach is more suitable 
for an infinite periodic structure. 

In general, the spatial periodicity is no longer 
preserved when the excitation is included in the formulation, 
especially when the excitation is a random field. In princi- 
ple, it is possible to express the total response under 
arbitrary excitation in terms of a fundamental solution which 
is the response due to one point-load on the structure. 
Mathematically, this fundamental solution is the Green's 
function and the total response can be represented as a convo- 
lution integral, but the actual calculation can become extreme- 
ly tedious. One type of random excitation which does not 
destroy the spatial periodicity of the system is that which 
is convected as a frozen-pattern at a given velocity 
Known as Taylor's hypothesis, this is an assumption freq.uentiy 
made in the analysis of i..irplane response to atmospheric 
turbulence. Unfortunately, significant decays in the corre- 
lation of pressure field have been found in experimental 
measurements of boundary- layer turbulences [^23-283, Thus, 
calculations based on frozen-pattern models are j'ust crude 
estimates as far as structural response to boundary-layer 
turbulence is concerned. Recently a scheme has been proposed 
C 29-303 in which a decaying turbulence field is constructed 
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from frozen-pattern components, thus retaining the computa- 
tional advantage of periodic structure models. This scheme 
(which we shall call the "turbulence decGmposit5.on scheme") 
will be used in this thesis and applications to structural 
response analyses will be discussed. 

More specifically, the advantages of the turbulence 
decomposition scheme will be discussed in Chapter II by 
comparing to the conventional point-load analysis. In Chapter 
III the basic concepts of the turbulence decomposition scheme 
and its application to the structural response spectrum 
calculation will be developed using a simple model of an 
infinite tmsupported beam exposed to a supersonic boundary- 
layer turbulence excitation, The effects of surrounding 
fluid will be ta^en into consideration. In Chapter IV the 
unsupported beam will be replaced by an infinite beam on 
evenly spaced supports which is a more realistic model of an 
actual aircraft panel system, The solutions will be compared 
with the experimental results. In Chapter V the case of a 
membrane with random structural properties will be considered 
and the effects on the spectral densities of the structural 
responses will be discussed. With potential applications to 
shin- friction drag reduction designs in mind, Chapter VI will 
be devoted to the analysis of the Reynolds stress in the 
boundary layer. 
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II. SPECTRAL ANALYSIS 


2 , 1 Introduction 

Measured cross-spectra of a turbtilence field usually 
show some decay in the statistical correlation in addition to 
convection at a characteristic velocity C23-28J, Under such 
a random excitation the computation of structural response 
statistics becomes much more tedious than that which would- be 
the case if the turbulence were coovected without decay; i.e,, 
convected as a frozen-pattern El^-193. The conventional method 
of analysis is a point-load approach. -.As it will be shown in 
the next section, this method requires a numerical double 
integration to compute the cross-spectral density of the 
response of a one -dimensional structure. If the forcing field 
is a convected froaen.-pattern field, then an alternative 
formulation will allow the cross-spectral density of the 
structural response to be computed v/ithout numerical inte^a- 
tions. However, because of the spatial decay in the measured 
turbulence spectra, the analysis for structural response based 
on the frozen-pattern assumption is just a crude estimate. 

The method to be discussed in this thesis retains the maximum 
computational benefits of the frozen-pattern analysis but at 
the same time the actually measured spectrum of turbulence 
field can be incorporated in the calculation. In this method 
a decaying turbulence will be treated as a superposition of 
frozen-patiTern rcomponents so that the structural response can 
be superposed simila'. ly. In the case of a one -dimensional 
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model, -this method only 'a sin^l^ inliegratiGn instead 

of a double integration in the point-load analysis, The method 
will he called the ’'turbulence decomposition scheme” in the 
sequel. 

To extend the turbulence decomposition scheme to a two- 
dimensional problem, the actual turbulence field is divided 
into strips which are running parallel to the x-axis, the 
direction of the mean flow. Each strip can then be decomposed 
into iiifinitely many frozen-pattern components as in the one- 
dimensional case. The total structural response is obtained 
by summing up the responses due to all the loading strips. 

In essence the point-load scheme is used in the across-flow 
direction whereas the frozen-pattern scheme is used in the 
along-flow direction, 

. In the following sections the one-dimensional case will 
be discussed first and then the two-dimensional structure. 
Furthermore, a general treatment of random structural properties 
will be included in the formulation of the two-dimensional 
problem, 

2,2 The Statistical Properties of the Excitation 

Measured frequency cross-spectra for pressure varia- 
tioi^ in a turbxilent boundary layer with respect to a fixed 
frame of reference have the general form of 


( 2 . 1 ) 
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where is a characteristic velocity, and and are non« 
negative definite even functions of % and , respectively* 
fhese two functions have an absolute niaximum eq.ual to one at 
the origin and they approach to zero at large absolute values 
of the argument, ^Che general form. Eg., (2,1), is sometimes 
attributed to Corcos £Z32» A number of researchers have 
reported curve-fitted results for $p{0,0,«^), and 

For representative works we cite the papers by Bull 
Willmarth and Wooldridge [25], and Maestrello, et al [26-283, 
Implicit in Eq., (2^1) is that a real turbulence is not a frozen 
one. We note that the above pressure spectrum reduces to that 
of a plane wave field if ^1^2 ('7 ) = and it would correspond to 
a frozen-pattern turbulence if “ i* 

For additional physical insight and later use in sample 
calculations, two measured spectra are given below* 

For a subsonic boundary'’ layer [26,273* 

■ - - 


_ 1 S* ^ ~ 


( 2 . 2 ) 


= exp C - ) 

0 

= exp ( - ) 

- 0 ,2if0 
A £ = 1,08 
= 1,80 
0^.-2 


(2.3) 

(2.4) 

= 0.470 

Kg 3.0 

K^ = 14, G 
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-3 -3 

0 = - d,2i|' X 10 ) + ltl5 X 10 sec 

V ^oc 

where U^o is the free-stream velocity, denotes the speed of 
sound in the field, S* represents the houndary-layer displace- 
ment thickness, and © is the average eddy lifetime. 

For a supersonic boundary layer £28] » 


— 1 s ^ 

® (0,0,W) 2 A 

P 2 n=l ^ 

- f- Knl“l 

(2.5) 

i)Tj( 6) - exp ( - ) 


(2.6) 

ilTgC’?) = e*P ( 

= h’tk' X 10 ^ 

-2 

= 5-78 X 10 

(2.7) 

-2 

Ag = 7.5 3C 10 

~ 2,k3 X 10 


-2 

= -9.3 X 10 

= 1,12 


Aij. = -2.5 X 10“^ 

= 11.57 


3 

0^2 = 0,26 



tig = 0,75 Uoo 

VIhere S is the boundary- layer thickness which is defined as 
the distance from the boundary at which the average value of 
the turbulence velocity reaches 0,99 Uoo* 


2,3 Poinx-Load Analysis 


Gonsider a unit concentrated load at x = $ on a one- 
dimensional structure where x is measured along the structure ^ 
parallel to the flow direction; We obtain a frequency response 
function to) by solving the equation 


JL C H(x,?,o>) e } = S{x-^) e 


( 2 , 8 ) 


where JIC } represents- a linear differential operator in x 
and t, and S( ) is a Dirac delta function, fhen^ the cross- 
spectral density of the structural response ®„(x^,x„,w) is 
calculated from 

ISL 

= jj H»(Xg,?g,lO) 




(2.9) 


Where the. cross-spectrum of the pressure field, 

i is the length of the structure under consideration, and an 
asterisk denotes the complex conjugate , It is assumed that 
the loading is a spatially homogeneous random process; there- 
fore, the cross-spectrum depends on Although this 

method is quite straightforward, long computer time is required 
to carry out the double integration in Eq, (2*9), 

2,4 Prozen- Pattern Analysis 

If the pressure is truly of a frozen type and is 
convected at a constant velocity in the positive x-direc- 
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tion, then it is a random function of x ~ Uj,t, Such a random 
function can he expressed as a Fourier-Stielt^es integral as 
follows s 


p(x 




i(f*3t - lex) 
e 

— Qd 


dF(k) 


( 2 , 10 ) 


wSere > the frequency and the wave-numher k are related to the 
convection speed as o>/k ~ It is known from the random 
process theory that 


S£.dFClCi } dF^ C^3 ) } - Sp(kj^ ) S C%-k 3 1 ^^2 \ ^ ^ ^ 

where E{ } represents the e^isemhle average, and Sp(k) is the 
wave-numher spectrum in a coordinate frame moving at the 
velocity (referred to as the moving frame in the sequel). 

The cross-correlation function E£p(Xj^-U^,t^) p(x 2 «IJQt 2 )} 
of the pressure, referred to the fixed frame, can he calculated 
simply hy use of Eqs, (2,10) and (2,11), This function, 
denoted hy Rp, depends only on S - U^T where ^ " ^2 

T = tj^ “ ’^2» related to the moving-frame wave-numher 

spectrum Sp(k) as follows 1 


EpC? - = 


f ik(Uo^ - ? ) , , 

e S„(k) dk 


( 2 . 12 ) 


* If a Riemann-Fourier transform is taken of Eq, (2,12) we 
obtain the fixed-frame frequency cross-spectrum of pi 
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®p(5.(o) 




Ep(§ - V^) 


= i_ s (ia.) 
Pci 


( 2 . 13 ) 


Equation (2.13) shows that the fixed-frame frequency cross- 
spectrum of a frozen-pattern turhulence has a special form 
■where ^ appears only in the imaginary exponent, Ehis equation 
also provides a simple formula to convert from S (k) to 

Jr*- 

$ Conversely, to convert from to S (k) ; 

P P P 


s^Oc) =|uj^^(o,ku^) 




Evaluated at ^ - o the cross- spectrum ^p(§,(f3) reduces, of course, 
to the usual spectrum. We' emphasize that Eqs,, (2.13) and (2,l4) 
are valid only if the turhulence is strictly of a frozen- 
pattern, and is convected at speed U^, 

Equation (2,10) suggests that the structural response 
to a frozen-pattern turhulence can he constructed from a 
fundamental solution where the excitation is just a convected 
sinusoidal pattern of unit amplitude, Thus, let H(x,k) exp(icot) 
he the steady-state solution for 


Ji{H(x,lc) exp(i(jct)} = expCi(wt - lex)] (2.15) 


Of course, this solution must satisfy all the necessary 
boundary conditions. Then the solution to 

{w(x,-t;)} = p(x - Ugt) 


(2.16) 
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after reaching stochastic stationality, may he expressed as 

OO 

w(x,t) = j H(x,k) exp(ioit) dF(3c) 

*«oo 

= I H(x,lc) expdUgkt) dF(fc) (2.1?) 

— OO 


It follows that the cross-correlation fiinction of the structural 
response is 


E£w{x^st^) wCxg.tg)} 
oe 

= J| H(x^,k^) e 


iU^(3Cj^ti — 


^ oo 


S S(fc^ - tg) dk^^dlCg 


OO * 

= ( H(xi,lc) H*(X2,k) SpOc) dfc (2,18) 


As expected, this correlation function is dependent only on 
t^ - tg. If it is desired to calculate this correlation func 
tion in the frequency domain, we may substitute Eq, (2 ,14) 
into^Eq.^ (2.18) and change Uj^k to to* 


E£w(x^,t^) wCXgftg)} 


( ia)(t- - to) 

J H(x^,£0/U^) H*(xg,w/U^) e ^ ^p(0,(d)du) (2,19) 

ecr 


In terms of the input and output spectra the relations are 
extremely simple and illuminative? they ares 
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in the wave-numher domain i 

S^(Xj,X2slc) - H(x^,k) H*(x2fk) Sp(k) (2.20) 

in the freq,uency domain i 

“ H(x^,co/u^) H*(x2>iaJ/U^j) ^p(o,a)) (2,21) 

VIhen XjL = X£ these formulas reduce to those for the usual 
spectra, and they have the same form as the well-known result 
for a single degree of freedom system in the random vibration 
t^heory. 

It is appropriate to call the H function in Eqs. 

(2,20) and (2,21) the wave-number response function for 
convected frozen load to distinguish it from the II fimction 
in Eq. (2,9) which is the frequency response function for 
point load, 

0?he advantage of the frozen-pattern assumption is 
clear, fo obtain the cross-spectrum of the response, no 
integration is required under this assumption while a double 
integration is needed in the point-load analysis, Eq, (2.9), 
Unfortunately, significant decays in the correlation have 
been found in experimental measurements of boundary-layer 
turbulences, Thus, correlations based on frozen-pattern 
models are just crude estimates, at least for the calculation 
of structural responses to the boundary-layer turbulence. 
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2,5 Turliulence Decomposition Scheme - One-Dimensional Case 

In the one-dimensional case measured frequency cross- 
spectra have the general form of 


= 0^(0, W) i}r(|) exp(-i<h^/u«) 


( 2 . 22 ) 


•which is reduced from Eq. (2.1), Implicit in Eq, (2,22) is 
that a real turbulence is not a frozen-pattern one, 

lo obtain a theoretical spectrum consistent with Eq, 
(2,22) the following representation of a general turbulence 
pressure is proposed* 


p(x 


,.)=J 


p(x - ut) dG(u) 


(2.23) 


Eq, (2,23) implies that p(x,t) is a superposition of infinitely 
many frozen-pattern components, each having a real random 
amplitude dG(u) and a eonvection velocity u, Suoh velocities 
can assume either positive or negative values. Of course, 
each frozen-pattern component can, again, be decomposed into 
sinusoids, ^IJhus 


p(x 


DO 

,«-JJ 


i(ufct - lex) 


dP(k,u) dC(u) 


(2.24) 


and its fixed-frame correlation is 


E{p(x^,tj^) pCxgttg)} 
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. CW 

= {[j] expCKuliiti 


Uglcgtg) - - ^2^2) 3 


ECdF(kj^,u^) dF*(i;;2»^2^ dGCug)) (2.25) 

In order that this correlation function may depend only on 
^ ~ ^ Xg and x: - t^ - tg I which we shall assume to he true, 

the ensemble average under the integral sign in Eq.. ( 2 * 25 ) 
must have the form 

E£dF(lCj^,u^) dF*(k 2 .u^) dG(u^) dGCu^)} 

= Sp(lc^,u^) - IC 2 ) S(u^ - U 2 ) dkj^d!3:2<i%<in2 (2,26) 

Substitution of Eq.. (2,26) into Eq., (2,25) results in 

Hj,(?,T) = j| SpC]c,u) di du ( 2 . 27 ) 

-*oo 

We now apply a Fourier transformation to obtain the fixed- 
frame frequency spectrum 

oo 

=2¥| e~^iX 

—OO 

= Sp(f,u) du (2.2B) 

Clearly Eq. (2,28) is a generalization of Eq. (2,13), 

To compare Eqs. (2,28) and (2,22), the latter is 
Fourier-transformed to yield 







OJ 


21T i *p 


U. 


(2.29) 


where 

?(v) 

Therefore , 


= 3f J >^(5) dg 


(2.30) 




$„{5.w) = $„(0,(J0) I T(o 4 - -ffi-) dJ. 


P 


(2.31) 


Letting ©4 = Ci/U| we obtain 


f 


®p(e.03) = *p(0,(U) 


/ 

— oo 


(2.32) 


U‘ 


IPheh -eouating and we find a foieiaula. to compute 

P P 


Sp(u>/u,u) as follows* 


= lul ^(f - f > 




= }g^| ip(0,(o) j <I>(S) exp|;i5(f - -^)3 d6 (2.33) 


— OO 


The freg,uency cross-spectrum for the structural 
response can be obtained by a similar superposition. Thus 
by a generalization of Eq.. (2.21), 


OO 

/ 




|u| 


H(x. 




H*(X2.^) 


SpC^.'^) 


du (2.3^) 


Or, letting t = 
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oo 


H*(x 2«^J Sp(fc,^|) di (2.35) 


Nowt since 


S_(lc,0 = |t| ! ® (O.to) If (kL -m 

P it P Uq 

we oTstaln a very simple result 

OQ 

:= lp(0,£i>) I HCx^.Jc) 

-^Oo 

K*(x2,k) ¥(lc .p: (2.36) 

As a check we note that when the turh^llence is frozen-patterns 
“ 1 and 

T(v) = S(v) (2.37) 

then SqL« (2,36) reduces to the same form as Eq., (2,21)* 

In the case when the structural acceleration response is 
the main concern of the problem, Sq,. (2,36) is changed to 

oa 

^ $p(0,U3) I H(x^,k) 

—00 

H*(X2,1!:) T(lc - ^) 6fc (2.38) 

Z^6 Turhulenca Decomposition' schema - Two-Bimensiorial- Case 

Generalization of the restilts obtained in the last 
section to a two-dimensional pressure field is straightforward. 


i? 


Assuming that such a pressure field can he decomposed into 
frozen-pattern components, we have 


p(x,y,t) 



ut,y) dG(u) 


(2.39) 


where both p and G are random functions and $ is a frozen- 
pattern component of p. Each component pressure ^ can again 
he constructed from frozen-pattern sinusoids, HJhus 


P(x 


oo 

.y.t) = jj 


i(ukt - hx) 


dF(k,u,y) dG(u) 




In order that the cross-correlation of the random pressure 
E{p(Xj^,y 3 ^,t^) P(3£^2»y2®^2^^ dependent only on ^ 

■J'y ss *“ ^2* ~ "^1 “ ’*'2* cross-correlation of dFdG 

must have the form 


E{dP(k3.,i;U^*y3^} 4F*(k2rW2»y2^ dG(u^) dG(u2)} 

=■ S(k^ - kg) S‘(uj^ * ^2^ 

dk^ dkg du^ dUg (2,41) 

Then, it can he shown that the cross-spectrum of p has the 
form 

^p(^i^»o)) J p S^(^,u,7j) du (2.42) 

>-oo 

Eq.uating Eq., (2,42) with the general form of the measured 
turbulent pressure spectra, Eq, (2.1), results in 
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Sp(^,u.i,) = - f-) tz(^) (2A3) 

o 

where is *fche Fourier transform of ^jr^i ige,, 

Oo 

?l{v) = 2W j d& (E.W 

— oo 

It is clear from Eq, (2,40) that, in this case, the 
structural response can also he constructed from a fundamental 
solution which is obtained hy letting the excitation he some 
suitable frozen-pattern sinusoid. Denote this fundamental 
solution by H(x,]c,y,y* ,co) and let H exp(ioit) be the steady- 
state solution for 

itot i(wt — fcx) 

X£H(x,3c,y,y*,a3) e } = e. S{y - y*) (2.^5) 

where JL represents a linear operator in x, y and t. She 
physical meaning of H is self-explanatory? it is the complex 
amplitude of the steady-state structural respor^e at (x,y) 
due to a strip of excitation, along y = y* which is convecting 
in the x-direction. Then, the solution to 

pC£w(x,y,t)} = p(x,y,t) \(2>/{-6) 

after 1 ‘eaching the stochastic stationality, may be expressed 


w(x,y,t) 



iojt 

H(xa3s:,y,y*,(A)) e 




as 
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b 

s f dy* 


oa 

( f 


ikut 


K(Xfl!i,y,y' ,lcu) e dFCk,u,y*} dG(u) (2. ^7) 


where b is the width of the structure in the y-direction, The 
cross-correlation of the response is, therefore, obtained as 
follows I 


E£w(x^,y^,tjL) w(x2,y2>*fc2)^ 
b. 


«*> 


ay| 




= J{ 

0 

® i i i 222 


Sp(lCi fUji^ ,y£ — S(h^ “ ^2^ S(u^ — U 2 ) 




j dyidy^Jj H^,k,y,,y£.fcu; H*(K2.k.y2.y|,lcu) 


ihut: 

e S (ls:rU,y* - yl) dk du 


{2.i^8) 


p' ' '-'1 '^2' 

Substituting Eg, (2,43) into Eg. (2.48) and applying a Fourier 
transformation, we obtain the cross— spectrum, of the structural 
response 


= ^p(0,0,tr) 


f 

dy£dy^ i H(Xj^,k,y^,y£,u)> 

J 

•^OO 
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H*(Xg,k,yg,y^,ai) ,W^(k - -^) fg(y^ - yp dk (2,^^9) 

She reason for not using the same decomposition scheme in the 
y-direction as that for the x-direction is the absence of true 
trend of convection which can be seen in Eq,. (2,1), Therefore, 
the correlation length of the random forcing field in the., y«r. 
direction seldom extends beyond one panel and the structural 
response can be computed q^uite accurately using a structural 
model consisting of just one row of panels running in the x- 
direction. For such a model the usual levy series represental. 
tion of the y-direction response is adequate and the double 
integration .on y£ y|- in Eq,’ '(2,^9) can. be. carried out 
without difficulty. 

Thus far, the structural properties are assumed to be 
deterministic; therefore, the frequency response function H 
is a deterministic function. Sometimes it may be of interest 
to include the effects of random variability of the structural 
properties in the calculation, then H becomes a random process. 
In this case the cross-correlation of the structural response 
becomes 


E£w(^,y3^,t^) w(x2,y2,t2)} 


)) 

G 




)) 


E{H(x^ * » y^ .y^ » ) 

>^2 *^2 dF(lc^ ,u^ ,yj^) dF'^^'ChgjU^jyj^) 


4S(«i) dG(u2)} ‘ WV 


(2.50) 


21 


If the structural properties and the turbulence properties are 
independent of each other, which is a reasonable assumption, 
the ensemble average inside the integral of Eq. (2.50) is 
separable. Again, using Eq. (2.^1) we have 


E{w(x^,y^,t^). wCxg.Ygitg)} 
b 

= 11 dy|dy^ J ECH(x^,fc,y2^,y^,lcu) H*(x^,fc,y2,y|»tu)} 

(2.51) 


, ikuz: 

Sp(k,u,y£ - y^) e dk du 


The cross-spectrum of the structural response is obtained by 
a Fourier transformation of Eq, (2*51)* resulting in 




= V° 


b 

,0,w) fj 

0 


oo 

/ 


dy:[dy^ 


E{H(x^,k,y^,y^,co) 


H*(x2,k,y2,y^,w)3 ^y^(h - ^) - yp ^ (2.52) 

2,7 Conclusion 


The turbulence decomposition scheme has been discussed 
and compared with the conventioimil point-load approach. It 
has been shown that a decaying and convecting turbulence 
pressure field can be constructed from frozen-pattern components, 
each having a different convection velocity, and that this new 
scheme simplifies greatly the analyses of random properties of 
structural responses under the excitation of boundary-layer 
turbulence pressure. In paticular, only a single integration 
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is req[uired to compute the cross-spectrum of the structural 
response instead of a douhle integration in the point-load 
approach, The integration involves a mve-numher response 
function which is the response of the structure to a unit 
convected sinusoid. In the following chapters applications 
of this scheme will he discussed for some structural configu- 
rations and flow field problems. 
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III. UNSUPPORajED INFINITE BEAM 

3.1 Introduction 

A simple model which, nevertheless, retains most 
import^t f ea.tures of aircraft panels in a houndary-layer 
environment is the infinite beam shovm in Pig,3<l* 
beam is hacked on the lower side hy a space of depth d which 
is filled with an initially quiescent fluid of density and 
sound speed a 2 » On the upper side the heam is exposed to the 
excitation of a supersonic houndary-layer turbulent pressure 
p, The fluid on the upper side of the beam which carries 
the turbulence has a free-stream velocity Up,, density and 
sound speed a^. Anjactual panel system of. an aircraft is 
reinforced by stringers and frames so that this unsupported 
beam is not a good representation at low frequency range. 
However, at high frequencies, the turbulence eddy size is 
much smaller than individual panels. Then the effect of the 
constraints at the stringers and frames becomes negligible,. 

In any case the infinite beam model is an ideal one which can 
be used to show the utility of the turbulence decomposition 
method as well as the structure -fluid interaction without the 
burden 'of mathematical complexities. 

3.2 Wave-Number Response Function of the Infinite Unsupported 
Beam 

As the beam responds to the excitation pressure p(x,t.) 
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its motion will generate additional pressures in the fluid 
media on the upper and lower sides. Denoting these generated 
pressures hy p^ and P 2 » respectivelyg the governing eq.uation 
of the heam motion is given by 


El 




^2 


° ^ - Pa)zi=o 


(3.1) 


Where E denotes the Young’s modulus, I is the moment of 
inertia, m is the mass per unit length of the beam. 

For the purpose of determining the wave -number 
response function, H(x,k), the turbulent pressure p should 
be replaced by expCiC^t - lac) 3 structural response w 

equated with H(x,lc) expCiwt) - A{h) exp(-ilcx) exp(ioit). 
Furthermore, we shall make the usual approximation that p^ 
can be calculated without regard to the presence of the 


turbulence. Then p^ is governed by the equation 




ai'^C 


2 ^2 
2r^ , ^ 


-) = 0 


(3.2) 


- 

and subject to the conditions that p^ can propagate only in 
the positive z-domain, and that 




AOc) 


i(cot - lex) 


'z=0 


'1 'st 


^x 


(3.3) 


The solution for p^, when evaluated at z = 0, is known to be 

Di] 

AOc) 


(3c4) 
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Some comments alsout Eq., (3A) are in order* (1) 0)/k is the 

speed at which the structxaral motion A(lc) expCi(ait - kx)] is 

propagated along the team, (2) A structural motion generates 

no pressure in the adjacent fluid medium if it is propagated 

at the same velocity as that of the fluid medium (the case of 

= Uw), (3) Theoretically, the generated pressure attains 

an infinite amplitude when the propagation velocity of the 

structural motion relative to the medium is eq.ual to the speed 

of sound (the. .case of jw/k - Uool = a^, the shock-wave effect), 

(i}-) When this relative velocity is less than the speed of 

sound! i,e., Iw/fe - Uo.,1 < a-* » the generated pressure should 

provide additional inertia for the structural motion (the 

apparent mass effect); therefore, a negative imaginary value 

2 2 i 

should he given to the square-root C(<n/k - Uoo) - in the 

calculation. 

The pressure generated on the lower side of the "beam 
is governed hy the equation 


3t^ 


2 , 3 ^ ■ , 

^2 2 ^ 2 *' ^*2 

3x 3z 


= 0 


and subject to the conditions 
3P2 


3z 


= 0 


at z ~ -d 


(3.5) 


(3.6) 


and 


3z 


a I? \ - lex) 

>(^ 2^0 A(k) e 


at z = 0 


(3.7) 


The solution for P 2 t when evaluated at z = 0, is given hy 
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, T ./>. cot(Yd) . itujt « lEx) 

<P2>z=0 * 5^ A(k) e 


(3.8) 


where 



(3.9) 


2 

For small d (shallow cavity) and Y >0, this pressure provides 

additional stiffness on the structural motion hut for certain 

ranges of d value it can change to an added mass, When 
2 

V < 0, Y hecomes imaginary in which case 


cot (yd) cothlYdl 

V " " Irl 


(3.10) 


Again, the term has the effect of an added mass on the 
structural response, 

Eq.s, (3.^) and (3.8) can now he substituted into Eg., 
(3.1) to find A(k) and, therefore, H(xjk), recalling that p 
must be replaced by expCi(wt - kx)] and w by A(lc)exp[|i(u3t-kx)3. 
3}he result may be expressed as 


H(x,k) = A(k) e 


-ikx 


r ^ 2 . 

= (Elk - mw + i^^a^ 




k(w/fc - U„)' 


C(ojA - 


•ikx 


(3.11) 


This wave-number response function H can now be substituted 
into Eg., (2,36) to obtain the cross-spectral density of the 
structural displacement response. 
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3.3 Numerical Example 

Numerical computations have been carried out for the 
freq,uency spectrum (i.e., when x^ = in Eq., (2*36) ) of the 
structural response using the following physical data* 
properties’ of ,the be ami 

El .(.bendirig rigidity) = 3.935 x N-m^ . 
m (^ss per unit length) = 9.7^5 Kg/m 
properties of the surroxindir^ fluid media* 

= f 2 = P (air density) = 0,11015 Kg/m^ 
ai = = a (speed of sound) = 261,6 m/sec 

lio (free-stream velocity on upper side of beam) 

- 575.6 m/sec 

d (cavity depth) = 0,1178 m 

properties of the supersonic boundary-layer turbulence 
pressure [28]* 

1 (0,10) - spectral density 

i(f(§) = decay factor = exp 

U^‘:( characteristic convection velocity of the 
ttarbulence ). = 0,75 

mfv ^ oj-) = — — Ir y - 

+ (h - 

S (boundary-layer thictoness) = 0.279 in 
experimentally determined constants 
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c< - 3 

4.2 

X 10 

Ag = 7.5 X lO"^ 

Ag = -9.3 X 10“^ 
-2 

Ajj, = -2.5 X 10 


■% = 5.78 X lo" 

K2 - 2,i|*3 X 10"^ 

Kg = 1,12 

Kij. = 11.57 


Fig, 3,2 shows the computed displacement freq,uency 
spectrtan of the structural response under the assumption of 
a truly frozen-pattern turbulence (i|y($) = 1), There appears 
only one peak around 475 Hz and the value of spectrum decreases 
rapidly as the frequency increases. The existence of the peak 
can he ea^Jlained as follows* In the absence of surrounding 
fluids j the wave number of the free structural motion would 

2 i 

be (moj/EI) , The wave number of the frozen-pattern turbulence, 
k, is related to the circular frequency oa by k = c^Uc. A 
•resonance, (called coincidence) occurs when these two wave 
numbers are equal £32], In the present case this coincidence 
frequency is found at 466 Hz, Therefore, the peak in Pig, -3.2 
results ftom the coincidence of the wave numbers of the frozen- 
pattern turbulence and the free structural motion. The small 
difference between the peak frequency in Pig, 3,2 and the 
estimated value above comes from the effect of the surrounding 
fluids on the structural response, 

Fig. 3,3 shows the spectrum of the structural displace- 
ment response when the measured spectrum of the turbulence 
pressure field is used in the computation. This spectral 
density has many peaks in contrast with only one peak in the 




frozen-paliterii case. Since infinitely many frozen-pattern 
components are present in the turbulence, the turbulence 
spectrum has also infinitely many mve-number components 
at each frequency. Thus, the coincidence (if we may use 
this terminology in this case) occurs at each frequency. 
However, the response magnitude may be small if a frozen-' 
pattern component which generates a coincidence wave in the 
structure has little contribution to the turbulence spectrum 
at that frequency. Therefore, even the coincidence is present 
at each frequency, its mgnitude varies as the frequency 
changes and some of them appear as peaks in the response 
spectrum, 

A comparison of Figs, and 3,3 shows that the • 
frozen-pattern assumption is unconservative which leads to 
lower estimates for the structural response and the radiated 
noise level 

3t^ Conclusion 

The theory developed in Chapter II has been applied 
to the simple escample of an unsupported beam exposed to 
boundary- layer excitations. The effect of a cavity and the 
effect of the free-stream velocity are included in the 
analysis. The spectral density of the structural displace- 
ment response was calculated usiaig the measured turbulence 
pressure spectrum of a supersonic boundary layer and the 
result was compared with the solution obtained for the ideal 
frozen-pattern turbulence. This comparison has shown that 
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the frozen-pattern assumption leads to considerably lower 
estimates for the structural response, Thus, we may conclude 
that the frozen-pattern assumption should not be used in the 
computation of the structural response spectrum under the 
excitation of a boundary- layer turbulence pressure. 



IV. PERIODICALLY SUPPORTED INPINIIE BEAT4 


ij- , 1 Introduction 

In this chapter the turhulence decomposition scheme 
will he applied to the analysis of an infinite . "beam supported 
at uniformly spaced intervals hy elastic springs, 5!he elastic 
supports are simplified versions of reinforcing stringers of 
an aircraft fuselage. Although an aircraft fuselage is a very 
complicated multi-panel system its dynamic behavior iss similar 
to that of the periodically supported heam described above. 

^Dhe one-dimensional beam problem, however, is more suitable 
for fundamental studies since basic concepts can be developed 
without the burden of mathematical details. Thus, the analysis 
of the present chapter will he restricted again to one spatial 
coordinate. 

At the first sight the problem of a periodic beam may 
not appear more difficult than that of any other structure if 
one accepts the linearity assumption and uses a normal-mode 
formulation. In practice, however, the normal mode of a 
periodic beam of many spans cannot be calculated accurately 
due to close clustering of natural freq.uencies in frequency 
bands. The futility of the normal mode approach in dealing 
with a large number of spans has led to two alternatives* the 
wave propagation approach (space -harmonic analysis) £17-211 
and the transfer matrix approach £10-1 61, The two alternative 
methods are closely related, however. The so-called free wave 
propagation constants in the first method are the natural 
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logarithms of the eigenvalues of the "basic transfar matrix in 
the second method £ 222 , ^Uhe computational simplicity in "both 
methods is obtained by utilizing the fact that the entire 
system is composed of identical sub-units in the formulation. 
The fundamental solution required for the construction of the 
total structural response is one corresponding to the excita- 
tion of a frozen- pattern sinusoid. To obtain this fundamental 
solution the formulation will follow Mead’s wave propagation 
method, but will take into account the effect of free-stream - 
velocity on the same side of the turbulence excitation and the 
effect of a cavity on the opposite side of the excitation. As 
a numerical example, the spectral density of the s-fecuctural 
response will be computed and the results will be compared with 
experimental measurements, 

^,2 Wave-Number Response Function of the Infinite Periodic 
Beam 

A sketch of the structural model is shown in Pig, ^.1, 
surrounded by the same acoustic environment as that assumed in 
Chapter III. 

The governing equation of the beam motion not directly 
ever an elastic support is given by 

3^w S w 

D-^ + m^= p+ (Pi - P2)z=_d 

oX ot 

where D denotes the bending rigidity and m is the mass per 
unit length of the beam. The additional pressure fields 
denoted by and p£ are generated in the fluid media on the 




37 


upper and lower sides, respeatively, due to the 'beam motion* 
For the purpose of determining the wave-numher. response 
function, HCXek), the turbulent pressure p should be replaced 
by exp[li(£ut - kx)3 and the structural response w by H{x,k) 
exp(iayb), Thus, Eq., (4.1) becomes 


(D 




,2 loot _ ^i{Mt-kx). 
mo) H) e - e + 


(Pi 


■ ^2^-^ 
z-a 


(4.2) 


The forcing function expCi(t*it - }ex) 1 gives every span the same 
excitation but with a phase-lag jHq = fcX from one span to the 
next. In this sense may be considered as the imposed phase- 
lag of the excitation. To satisfy the spatial periodicity in 
the structural response Mead suggested the following series 
form Cl83* 

H(x,h) = S „exp(- - ^S ^ ) (4,3) 

n=-oo ** X 

where 

~ J^o (4,4) 

Without the elastic supports the wave-number response function 

would be just the one term associated with the forcing phase- 

lag yC/Q, The elastic supports give rise to multiple reflec-’ 

tions, thereby admitting other values. 

For the induced additional pressure p , we apply the 

1 

same assumption used in the preceding chapter. Then p^ is 
governed by the eq.uation 
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„ a .2 „ 2 .a^ ^ a ^ 

^ ^ ax^ % " ^-1 ^1 ^ 


(i^.5) 


az^ 


and subject to the conditions that p^ can propagate only in 
the region s > d, and that 


ap 




itot 


(^.6) 


^The solution for p^ corresponding to each component in the 
wave-number response function exp(-i^^x/-t) is known. The 
total can then be obtained by superposition, This pressure, 
when evaluated at a - d, is given by C31Ii 


Pllz=d = ni-«^ 

i((ot - /»n^ ^ 


(u^ - UJ^ 

W/5Un)C<“n - 




where 


wi 




The same comments which were made previously relative to Eq, 
(3.^) apply in this case to each term in Eq. (^.7), Each 
component now has a different propagation speed u^^ instead 
of the one propagation speed oj/k in Eq. (3»^)t emphasis 

and clarity these comments are x^estated as follows x (1) 
the speed at which the component structural motion expCi(o3t 
- is propagated along the beam* (2) A component 

structural motion generates no pressure in the adjacent fluid 
medium if it is propagated at the same velocity as that of the 


fluid medium (the case of ui^ = U^). (3) Theoretically, the 

generated pressure attains an infinite amplitude when the 

propagation velocity of the structural motion relative to 

the medium is equal to the speed of sound (the case of (u^^^ - 

Ucol “ “fche shock wave effect), (^) ’When this relative 

velocity is less than the speed of soundi i.e., ju^j^ - Uool < 

the generated pressure should provide additional inertia for 

the structural motion (the apparent mass effect) ; therefore, 

a negative imaginary value should he given to the square-root 
2 :2 ^ 

l-<"n - Ueo) - ] in the calculation. 

The induced pressure Pg in the fluid medium 0 :£ a < d 
is governed hy 


^2 - 

^ P2 ^ 2/a' 
2 > “^^2 VT 




hir 9x Bz‘ 

and subject to the conditions 




(^.9) 


= 0’ at a ” 0 (^#10) 

H at a = d (4,11)- 


The solution for pg, when evaluated at a = d, is given by 



ea 

s 

n=-eo 


COt(yytd) 


n 






(4.12) 


where 
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2 

For a small d and a positive , the n-th component of this 

transmitted pressure Pg gives rise to additional stiffness 

on the structural motion and for certain ranges of d value 

it can "become an added mass to the system. When T is 

n 

negative, "becomes imaginary in which case 


cotCr^^d) 


n 


coth|y^d| 



(4.14) 


and the component always contributes to the system inertia 
regardless of the value of d. 

Equations (4.3), (4.?) and (4,12) can now be substi- 
tuted into Eq, (4,2) to obtain 


oo 

^AjjCp(n) exp(-ijAj^x/£) = exp(-ijXQx/&) (4.15) 


■ where 


tp(n) = D (^) - 


(% - geo) 




2 oot(ynd) 


(4.16) 


n 


To determine the amplitude we follow Mead's 
procedure Ci9i] and calculate the rlrtual work done by the 
external forces acting on the structure and by the internal 
forces in the structure through a virtual displacement 



^-i(o)t - 


Excluding the elastic supports, the virtual work done within 
one span of the beam is 


X 





0 


(^. 17 ) 


The elastic supports are characterized by a translational 
spring constant K-^, a translational inertia M, a torsional 
spring constant K^., and a torsional inertia I, Thus the 
virtual work contributed by each elastic spring is 


SC 


2 


(K - SA- 2 + (kj. 

'n=-oo 


Io>^) 






Since the structural motion is spatially periodic the virtual 
work done throughout the entire structure is proportional to 
that of a periodic unit. Therefore, the principle of virtual 
work can be stated for a periodic unit as follows s 

+ SW^ + SW^ = 0 (^.19) 

which leads to the simultaneous algebraic eq.uations: 

2 2 
, ^ Kj. - Mw CO Kr - Ito 03 

A.cp(5) + ^ A„ H. — ^ ^ 
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(4p20) 


/ 1» if 3 - 0 
I 0, if j # 0 


In aci 3 ual compuijations the number of simultaneous equations 

must he truncated. One can, for example, solve a system of 

2N + 1 equations corresponding to -N ^ 3 N, The choice of 

N must he such that the truncated version of the wave-number 

response, Eq. (4,3), do not change appreciably by further 

increase of the number of terms used in the computation, 

2 

Equations (4,20) are derived for finite - Mw and 
finite Kp - lo?, and these equations cannot be reduced to 
those for supports rigid in translation or rotation. For 
example, if the supports are rigid in translation then 
becomes infinitely large, but the summation of all the 
must be zero since the deflection at each support is zero, 
and the product of these two becomes indefinite. To deal 
with this case, substitute 

An = - S A„ (4.21) 

n^^O 

into Eq. (4,3) to obtain 


H(x.k) = S A 
n=-«o ^ 


(4.22) 


Correspondingly, virtual displacements are chosen in the 
form of 

iM4x/!i -icat 

SA. (e - e ) g 


Then, instead of Eq., (ii-,20), one obtains from a similar 
derivation 


Aj<p(3) + cpC0)J_^Aj, 
n^O 

= - 1 , 3 # 0 


Ky - 103 M 

X X n=-^ I ^ 

n^feO 


(^.23 


It is interesting to note that if the supports are rigid in 
translation hut without constraints in rotation (the case of 
hinge supports) the equations for Aj can he decoupled* For 
such a case, Eq, (^. 23 ) reduces to 


AjCp(j) + = -1, 3 5 ^ 0 

Eq, (^,2^) shows that the product A.q 3 (j) is independent of 3 

tl 

e * , 

A^epd) = A2'P(2) = ..... (4.25 

Ihus, substituting 

= Aj(p(3)/q3(n) (4.26 

into Eq, (4,2^), one obtains equations involving only one 
unknown! 


A.:Cp(a)q3(0) S 

^ n--oo 


1 

cp(n) 


(4.27 


which is solved readily to give 


(4.28) 
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Uncoupled solutions such as Eq., (4.28) are not 

restricted to the case of hinge supports. In fact, if one 

of the two infinite sums in Eq, (4.20) can he dropped, a 

substitution of the type of Eq, (4,26) is possible which is 

the key for reducing Sq. (4.24) to Eq. (4,27) containing only 

one Aj in each equation, !This is the case when either 
2 2 

- Mw = 0 or - 1(0 = Oj i.e,, when the elastic supports 
offer no translational constraint or no rotational constraint. 
If the rotational constraint is infinite, one again 
cannot obtain a reduced equation from Eq. (4.20) which is 
valid only for finite support constraints. To derive a reduced 
equation one must use the zero slope condition at the supports 





(4,29) 


The remaining procedure is very similar to that leading to 
Eq. (4.24) 

Now, we are ready to compute the spectral density of 
the structural response. Once Aj are determined by either 
one of Eqs. (4,20), (4.23) or (4,28) depending on the problem, 
they can be substituted into Eq. (4,3) 'to obtain the wave- 
number response function, H. ' Further substitution of this 
H function into Eq, (2.36) along with the statistical proper- 
ties of the turbulence pressure gives the cross-spectral 


density of the structural displacement response d 
4,3 Numerical Example 

To illustrate the application of the present theory, 
the spectral densities of the acceleration response at a mid- 
span location (i.e,, - Xg = have been computed "based 

on the following physical data* 
properties of the heam* 

D ^{‘bending rigidity). ,= 3*935 x 10 N-m 
X (span-length) = 0,508 m 
m (mass per unit length) = 9*7^6 Kg/m 
properties of the. surrounding fluid media* 

Pi ” f 2 = P (density) - 0,11015 

a^ = a2 = a (speed of sound) = 26l,6 m/sec 

Uoo (free-stream velocity on upper side of "beam) 

= 575# 6 m/sec 

d (cavity depth) = 0,1178 m 
properties of the turbulent pressure C28]]i 

— ' 1 S ^ 

$ (0,o>) = spectral density = 'o 77" S 
P ^ Woo n=l ^ 

ij;(§) = decay factor = exp(--^) 

O 

(characteristic convection velocity of the 
turbulence) = 0,75 

S (boundary-layer thickness) = 0,279 ^ 



A-6 

-2 

A- == kA X 10 

-2 

Ag - 7,5 X 10 

A^ = -9.3 X lo"^ 

Ai|. = -2.5 X 10"^ 

Except for the additional information about the span length X 
the above physical data are the same as those used previously 
for the unsupported beam, and they were taken from a recent 
experiment on a multi-panel system £283. ^ke structural 
specimen used in this experiment was actually a two-dimensional 
panel array as shown in Pig. 4.2, Therefore, some data have 
been converted to their one -dimensional equivalents. For 
example, the bending rigidity D of the beam was the average 
value for the skin and the reinforcing stringers over a unit 
width, and the specific mass m was obtained similarly. 

Although the actual structural specimen had only seven spans 
and the two end-spans were somewhat shorter, it was felt that 
the theory of an infinite periodic beam on evenly spaced 
supports shotild give a reasonable result for the response 
spectrum at the center of the middle span where accelerometer 
A20 was located (referred to Pig, 4,2), and where the effects 
of the end spans were least important. The translational 
constraints provided by the supporbing frames were sufficiently 
strong to justify taking the translational spring constant ICj. 
of the supports to be infinite (i«e,, the deflections at the 
supports were assumed to.be zero). For the rotational con- 
straints we selected Kj, == 60 K-m/rad and I = 3.3 x 10“ Kg-m*^ , 
These are the one-dimensional equivalents of the torsional 


- 5.78 X 10 
K2 = 2.43 X 10 
= 1.12 
= 11.57 


"1 

yi 
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constraints of the frames if the torsional mode of each frame 
is a half-sine curve. 

In order to minimize the computer costs the computation 
was carried out only to 3000 Hz and at the intervals of every 
50 Hz, The resolution of the computed spectrum was compromised 
somewhat by the use of coarse intervals, but our main objective 
was to find the general trend which could be revealed by the 
values at $0 Hz intervals. 

In Fig, ^.3 the computed spectrum is shown along with 
the experimental spectrum. As it is customary, the experimen- 
tal speclmm is one-sided (restricted to the positive frequen- 
cy domain); therefore, the theoretical, spectrum has been 
converted by multiplying the computed two-sided values by 
two. It also should be noted that the experimental results 
were obtained using a filter of 1 Hz bandwidth. This accounts 
for its much more rugged appearance than the theoretical one 
computed at much larger intervals of 50 Hz. Furthermore, 
experimentally obtained signals may contain noise other than 
the structural response, ' Although the theoretical and the 
experimental curves show the same general trend, the former 
is lower than the latter throughout the entire frequency 
range investigated. This is to be expected since the theo- 
retical curve represents the average between the panel response 
and the stringer response, whereas the experimental curve 
shows the panel response alone. 



theoretical 


experimental 



frequency, Hz 


Fig, 4,3 Acceleration spectra of structural response 
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^.3 Conclusion 

The case of fluid- loaded infinite periodic beam, 
considered in this chapter, is one of the very few where 
a mathematically exact solution for the wave-number response 
function, H, can be obtained. If the beam is finite in length 
then the method of transfer matrix may be more preferable,' 
however, the effect of fluid loading cannot be accounted for 
exactly (in the mathematical sense) at the present time. 

Further extensions to the two-dimensional case of 
panel systems are obvious. If only one row of panels is 
considered, and if the two parallel edges of the panel row 
are assumed to be simply supported, then separation of spatial 
variables is rjossible in expressing the structural response. 
This is the well-knovm Levy's type solution for plate 
problems. With small modifications, the solution for the 
one -dimensional beam case can be changed to suit such a panel 
row problem. When more than one row of panels are included 
in the structural model the separation of spatial variables 
in the structural motion is no longer mathematically exact, 
but a separable form can still be used as an approximation. 
Although new concepts are not rei^uired in-:' treating such two- 
dimensional problems, the machine computation time can become 
extremely excessive and bm-densome to small research budgets. 
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V, MEMBRANE VaOIH RANDOM TENSION 
5 • 1 Intr oduct i on 

The unsuppor*bed infinite 'beam and the periodically 
supported infinite heam discussed in the preceding chapters 
are idealised models * However, actual sti'^uctures cannot he 
constructed in ideal manners. The material properties vary 
randomly throughout the entire structure and manufactiiring 
errors always exist in sise and shape such as the span length 
and the cross-section of a beam, the pre-tension in a membrane, 
etc. It has been shown by Lin and Yang X33- 351 that in the 
case of a periodic beam the randomness in the structure 
properties causes appreciable variations in the structural 
response from the ideal model results. In this chapter, the 
structural response of a membrane to the subsonic boundary- 
layer turbulence will be investigated. The pre-tension in 
the membrane v/ill be treated as a random process in space, 
and the membrane is surrounded by an acoustic environment 
similar to that in the preceding chapters, 

5,2 Wave-Niamber Response Function of Membrane 

The structural model chosen for the present study is 
a membrane which is infinitely long in the x-direction and is 
fixed along y « 0 and y = b as shown in Fig, 5,1. This 
membrane is backed on the lower side by a cavity of depth 
d which is filled with an initially quiescent fluid of density 
p 2 and. sound speed a^. On the upper side the membrane is 



Pig* 5*1 An infinite membrane with random tension under the 
excitation of boundary- layer turbulence 
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exposed to the excitation of a subsonic boundary- layer turbu- 
lent pressure p. The fluid on the upper side of the membrane 
which carries the turbulence has a free-stream velocity Uoo, 
density and sound speed The design pre-tension in the 

membrane is Tqj however, it will be assumed that this uniform 
tension can only be achieved in the x-direction but in the 
manufacturing process a small random variation €f (x) has 
developed in the y-direction where € « 1, 

As the membrane responds to the excitation its motion 
will generate additional pressures in the fluid media on the 
upper and the lower sides. Denoting such pressures by and 
P 2 » respectively, the governing equation of the membrane is 
given by 


m 




a w 


S W 


= P + (Pi - P2)„_, 


(5.1) 


where m is the mass of the membrane per unit area, Y is the 
viscous damping coefficient. 

As discussed in Chapter II, the turbulent pressure 
p should be replaced by expCi(wt - lcx)3 S(y - y*) and the 
structural response w by H(x,h,y,y* ,io) exp(itiit) for the 
purpose of determining the wave-number response function 
H(x,h,y,y*,ti^). Thus, Eq. (5.1) becomes 

2 2 

t iY( 0 H - Tq - [Tq + €f (x)3 

ax^ a^r 
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i{03t - kx) ,v ^ / X 

= e S(y - y*) + (p^ - 


( 5 . 2 ) 


The boundary conditions for the wave-number response function 
H are 


H(x,fc,y,y*,co) =0 at y = 0, b (5.3) 


These boundary conditions are automatically satisfied if H is 
expressed in a levy type series 


oo 

K(x,k,y,y',(b) sin j Aj^(^,k,y' ,£0) e"^^ djj 

(5.^) 

In the sequel the symbol A^^C^a) will be used in lieu of 
ICiy',to) for compactness. 

For the determination of the radiation pressure p^, 
we use the usual approximation that it can be calculated 
without regard to the presence of the turbulence. Under this 
condition* the radiated pressure p^ is governed by 


^Bt 


B 

^Bx 


) P. 





0 


( 5 . 5 ) 


with the boundary condition 



z=-d 



IT 

H e 


( 5 . 6 ) 


The exact solution for p^ is unknown but in close proximity 

to the membrane p can be approximated by the expression 

1 
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iwt “ . 

e S sin — — 
n-1 b 




(5.7) 


/ 

— oo 


Substituting E^i (5*7) into Eg.. (5*5) » multiplying both sides 
. 0 ^ the equation with sin(n* 7 ry/b) exp(iyji*x), and integrating 
over < X < CO and 0 < y < b| one obtains 

^(z,^) - fej3^(/i) gjjCz,^) = 0 (5.8) 

where 

=/ + ( 5^)2 . ( 5 . 9 ) 


Equation (5*B) is solved, with the condition that only 
propagates in the negative z-domain, to give 


g^(z,yi) ,= Cj^(/x) expClCj^(^)z] (5.10) 

Substituting Eqs. (5.^)* (5.7) and (5.10) into Eq. (5.6), and 
using the orthogonality properties, we have 




{fx> - 

!%(/) 






(5.11) 


Therefore, the radiation pressure p^, when evaluated at 
Z -d, is 


z=-d 


0 icot “ . ^y 

Ti e S sin —r— 

1 21=1 ^ 


' (lu - 






djx 


(5.12) 
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The transmitted pressure 
Helmholtz equation 


is governed by the 


.2 

^ P2 


a 

®2 I 2 

dX 





0 


and is subject to the boundary conditions 
^P2 

T — = 0 at z = 0 
3z 


^pg 

dz 



H e 


iwt 


at z 


-d 


SP2 

— = 0 a-t y = 0, b 


(5.13) 


(5.14) 

(5.15) 

(5.16) 


fhe conditions (5.14) and (5.16) corresponds to the ideal 
case where the bottom and side walls of the cavity are 
acoustically hard. The solution of this system is given by 


v/here 



i(iit fio « 

e S S cos 

n=l r=0 

n'Hr=odd 


rtry Zf n 

^ TTg(r) ^2 _ p2 


OO 

( cot lc^(/^)d 

J ^(/) 


V/*) 



(5.17) 



= 

(5.18) 


|2, r = 0 


g(r) = 

1 1, r ^ 0 

(5.19) 


Substituting Eqs, (5*^)j (5.12) and (5.1?) Into Eq. (5.2), we 
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have 

OO 

Ml 

OO ^ 

+ eQ^ j A^w c^ j f(x) e^^>* ■ dx] dX 

— “O Jtkei 

= (|-) S(^ - kj sin ^ (5.20) 

where 

cpj^(/i) = -mu? + iro) + + ‘^qI^ - fi 

(5.21) 

, , 2 ^ 16 n i 

"'nJl'^^ =<2“ rio ,r®s(r) n2 - r2 - r2 

n+r-odd 

i+r=odd (5.22) 

In the present study Eq,. (5.20) will be solved for the 
following two cases* 

case I I f(x) = 0 

U 

case II t f(x) = 11 X. exp(iaux), ^0 = 0 

d=-M 

Case I corresponds to the ideal uniform tension problem* i.e,, 
the tension in the membrane is eq.ual to the constant Tq in all 
directions. Case II represents the case where random varia- 
tion of tension in the membrane can be expressed as a super- 
position of sinusoidal variations, The fundamental wave 
number v will be assumed to be deterministic and the complex 
amplitudes I X::» random variables. 


(m - 
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CASE I I fUy - 0 


In ■this case Eg.. (5.20) reduces to 


oo 


niry' 


(5.Z3) 

which is an infinite set of simultaneous equations in A^. In 
actual computations the number of simultaneous equations must 
he truncated. One can solve a system of N equations corre- 
sponding to 1 ss n ^ N. The choice of N must he such that 
the truncated version of the structural response, (5.^)* is 
sufficiently accurate and it does not change appreciably by 
further increase of the number of terms. After -truncating 
the number of simultaneous equations, Eq* (5.23) can be 
written in a matrix form 


CKj^()j)] £a(^)} = Cp} S(^ - k) 


(S.z'f) 


where the elements of the matrix and the vector P are 




(5.25) 


^ .2. . W 

I’ = (tr) — r — 

n b b 


( 5 . 26 ) 


and is a IQponaolcer delia having tihe propertyt 


Snt= I 


1, n = i 

0, n ^ i. 


(5.27) 


The solution of Eq. (5.24) is 
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£A(j*)} = £p 3 - fc) (5.28) 

SulDstitution of Eq., (5.28) into Eq, (5«^) yields 

N uTty 

H(x,Ic,y,y',a3) == E R sin-— (5.29) 

n=i ^ ^ 


where is the n-th element of the vector H which is defined 
as 

£R} = £%(!:)] £p} (5.30) 

CASE II s f(x) = E exp(ijvx), Xq = 0 
In this case Eq, (5.20) reduces to 
sp^if^) A^(jA) 


9 nn V ‘ 

S X.A_(U + Jv) = (— ) sin-r— S(u - k) (5.31) 

j=-M D J 

where the numher of the simultaneous equations has "been 

truncated to N, In a matrix form 

M 

1;Ki(/*)3 £a{;^) 3 + eCKgJ .s ^ £a(ji( + -jv)3 


= £P} S(;i - k) 


where the elements of the matrix Kg are 

,mr.2 


(5.32) 


(5.33) 


6G -1 

'A 

1 

which show that the matrix Kg is a diagonal matrix, solve f| 

Eg, (5*32) £A^)} is expanded in a series, keeping in mind j 

'..j, 

that € is a small q^iantity, 

Ca^) 3 = + et^A^)} + 6 ^£^a^) 3 + ■ • • (s.S'f) | 

Substituting Eq, (5«3^) into Eq, (5*32) » and grouping terms 
of the first and second powers of €, one obtains 

CKi^)3 £°A(jh)3 = £p 3 - k) (5.35) I 

M '■ 

Dci^)] £^A ^)3 = - CKg] .s Xj£ A^ + jv)3 (5.36) I 

^ i 

2^1 1 

Ta^)} = - ^ A^ + jv)} (5.3?) 

Eq. (5.35) is solved to give 

£°a9j) 3 = C%9i)3‘^ £p 3 - k) (5.38) j 

Substituting Eq, (5«38) into Eq, (5-36) f we have 

V? 

, .. M 1 I 

£^A^;i) 3 = - D!i^)3 CKg] S ^ Xj[Ki^ + jv)]" 

•'T 

£p 3 + jv - k) (5.39) 

Likewise , 

1 M M i 

C A^)} = CK.(M)r^CK^3] S E X. Xei* 

^ ^ ^ 3=-M q=-M ^ ^ 

- qv)]'^CK2^ + jv - qv)] ^ j 


6l 


Cp} S(j 4. + - q.v - k) (5.^0) 

2 

Thus, to the order of € , we have 

^ iMiy ^ j 

H(x,k,y,yS(jo) = Z sin-r— CR j^ + € 2 ;P“ 

n=l ° j=-M ^ ^ 

2 ^ * 3<1 

+ e 2 2 X. X« G ] (5,/}.i) 

5=:-M q=-M 3 ^ n 

where has been defined in Eq., (5i30)f an<3L and are 
the n-th elements of the following two vectors, respectively, 

{P^} = -C%0c - 3V)f^CK2: CKia£)f^£P} 

CG^'^} = CKj^Oc - 3M + Gv):‘^£K 2D CKj^Ck - 3 v)3"^CK23 

(5.43) 

Now, if the random variables Xj, j =1,2,* • *,M, 
have zero means; i.e., 

eEx,] = 0, 3 =‘-M, • ■ 1, • • ■, H (5.44) 

V 

then the mean value of the wave-number response function H is 


eCh] = 


K 

2 

n=l 


sin 


mry 

b 




M M 
3«-M ci=-M 


ui 


Cx 


0 04 ^^ 


(5.^5) 

We note that the first term on the right side of Eq, (5*^5) 
is the same as the wave-number response function for the ideal 


uniform tension membrane, (5.29) » and the second term gives 
the contr?*bution from the random variation in the tension 
field, fhe cross*-spectrum of the response function H is 


BLH (Xj^ . h* p y^* , oj) H* , k, y^ , , 03 ) ] 

N N nTTy ^ 

“ S £ sin — — sin ^ — {R R« + € S 2 

n=l i=l ” t ^ n 1 


ECXj \1 C4 4* + E„ af * + (5A6) 


where and y£ should he substituted into the variables with 
the subscript n and and yj should be substituted into those 
with the subscript 1, 

For the special case where are statistically in- 
dependent of each other and identically distributed random 
var i able s , and 


ECXj X*] 


2 . _ 
r CT • 3 “ *1 

^0 * d ^ 


(5.^7) 


Eq.s, (5.^5) and. (5.^6) reduce to, respectively, 


N 


E[H] = S sin ^ Ci? ^ 
n=l ° ^ 


M 

MZ r. l ^ JjO 2 am 
2M j=-M ^ 


(5.48) 


N N rsjry- Jiry2 „ 

E[H H*] = 2 2 sin sin 

n=ll=l ^ b ^ X 
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where e®cr 


(5.50) 


is the coefficient of variation of the random tension in the 
y-direetioiio We note that due to the assumption of uncorre- 
lated X-, Eq., (5*^9) “becomes linear in each harmonic component 

u 

of random tension field in the sense that the solution can be 
obtained by solving for each harmonic tension field variation 
separately and then superimposed to obtain the total correla- 
tion, 

Nowi the cross-spectrum of the structural response 
can be obtained by substituting into Eq., (2,52)* Generally 
the integration over k in Eq« (2,52) must be performed numer- 
ically, but the integrations over y£ and y^ can often be 
calculated analytically. To illustrate we shall assume that 
the decay factor of the pressure field in the y-direction, 
^2^^i “ ^2^* expressed as 




( 5 - 51 ) 


vrhere Q is an experimentally determined turbulence scale. 
Substitution of Eqs, (5-^9) and {5* 51) into Eq, (2,52) 
results in 


V^’='2»yi>y2'“^ 


r 


== $p(0,0,co) 


N N rary- HttYo 

^ (}c - ^) £ S 2 sin = sin 

^ U« n=l 1=1 b b 


N W ” ^ 

S 2 []U_«(X4) oisji ^ 


6h 


■*• V,j^5 ^st (5.52) 


where U„_, V__ *, and W„_ :• are the elements in the n-th 
ns ns t J ns | j 

row and s-th coliimn of the following matrices, respectively. 


CTI(x): = (-I) CKi(k) 3 '^ 


(5.53) 


1 2 ' r*— y - • V .n**l — -i .» * . » -i"l 


CT.(x)3 = - (^) CKjdc - jv)3 "GKg] CKjCk)] 


J "b 

-i{k - jv)x 


(5.5^) 


i;wj{x)3 = (-I) [K^(k)r^[K2: ck^dc - 
CK^(k) 3 -^ e"““ 


(5.55) 


and 


I = ^ JL c 

(l/Q)^ + (sx/b )2 Q 

^ B 7 l/h tT/h 

( 1 /q)^ + (sjiA)^ (i/s)^ + (-tVti)^ 

£1 - C(-D® + (-1)"^] (5.56) 

In carrying out the numerical computations the most time- 
consuming part is the inversion of matrix heside the 
numerical integration over h which cannot be avoided, 
Therefore, in the actual computation, only the diagonal 
terms of matrix will be kept to save the computer time. 
Hence the coupling of modes due to the effect of transmitted 
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pressure, P2, will "be negleciied, 

5# 3 Numerical Example 

The following ph^^sicaX data are used in the numerical 
calculation; 

properties of the membrane (Mylar membrane )i 
p (density) = 1,304 x 10“^ Ib'-secVin^ 
t (thictaiess) = 0,00035 in 
m (mass per unit area) = pt 
To (tension) 0,35 lb /in 
b (width) = 7 iii 

Yqi (coefficient of the random tension variation) = 0,1 
V (fundamental wave number of random tension harmonics) 
= 7T in“^ 

properties of surrounding fluid media* 

pj = P2 “ density *= 1,15 x lO”*^ Ib-secVii^^ 

4 

“ ^’2 ~ speed of sound = 1,3^ 2C 10 in/sec 
U^(free stream velocity on the upper side of the 
membrance) = 1200 in/sec 
d (cavity depth) - 1 in 

* 

properties of turbulences c 

- , 1 S* ^ 

(0,0, co) = spectral density = "o 17“ ^ ’^n ® 

U (characteristic convection velocity of the * 
o 

turbulence) = 0,8 Uao - 1^1 

= decay factor in th^ x-direction = e 
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1 -n-tJ^jeCCl/Doe)^ + (k - 

ilr^CTj) - decay factor in the y-direction = 



(Q ** oi$*) 


S* (boundary- layer displacement thickness) 

= 0,0ij'6*x-Re"^*^ 

Re = l^x/v 
X = in 

V = 2,269 X 10 "^ in^/sec 
and experimentary determined constants 
= 0,240 = 0.470 

Ag = 1,08 Kg = 3.0 

= 1.80 = 14.0 

C3( = 2 

6 (eddy lifetime) = -(1,24 x 10“^} (Uoo/a^ ) 

-3 

+ 1,15 X 10 sec 


Fig. 5.2 shows the spectral densities of structural 
response at x^ = Xg (the spectral densities of structural 
response is homogeneous in the x-direction, hence it is not 
a function of x^ or Xg) and ® with or without the 

random variation in the membrane tension computed under the 
assumption of a frozen-pattern turbulence (i.e,, =1). 

Purthearmore , only two harmonic terms of the random tension 

—1 “1 

field, corresponding tov^irin ( 3 -I) and -V = -TT in 
(5 == -l)c have been included in the computation. We recall 
that the cross-correlation of the structural response is 
linear in each harmonic component of the random tension field 
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when the random amplitudes of the tension, X^, are statist!- 
cally uncorrelated with each other* The spectral density of 
the structural response in the absence of random tension 
shows no peak within the frequency range shown in the figure, 
while several peaks appear if the random tension is taken into 
account. This change in appearance cam be explained as 
follows. Since the exciting pressure is assumed to be a 
frozen-pattern turbulence, the wave number k of the turbulence 
is related to the frequency by k = Without random 

tension the wave number of the response,^, of an infinite 
membrane unsupported in the x-direction must be equal to 03/u , 
Therefore, the coincidence resonance C 32 D would occiar when one 
of the diagonal elements of the influence coefficient matrix 
could become zero. If the induced pressures, p^ and p^, 
are neglected, this would happen when 



However, numerical computations have shown that the left-hand 
side of this equation is always positive for any value of tu , 
Thus, coincidence can never occur when uniform tension Tq is 
acting on the membrane alone. When the random tension terms 
are included in the computation, there appear perturbation 
terms which have factors of the form of - 3v)j 9 j . 

t!,. » •, in the wave-number response function H, Thus 

a shift of the wave number in the structural response from 
k to k - occurs. When 3 = 1 , the coincidence frequency 
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can be calculated from the equation 

2 , _ ^nTT* 2 . m / ti) \2 _ 

- mo> + T^(~) t ^o^ir “ ^ 

c 

For n - 1 and 2, this equation has solutions 

f = 371 Hz, 720 Hz, for n = 1 

f = 425 Hz, 665 Hz, for n = . 

No solution exists if n is greater than two. There is no 
coincidence resonance when j = -1, Note that these values 
are only rough estimates of the peak frequency since the 
induced pressure fields, and p , which have been ignored 
in the estimates provide additional inertia or stiffness to 
change the structural response, thus altering these peak 
frequencies, 

Pig* 5«3 shows the spectral densities of structural 
response using the measured spectrum of turbulence pressure 
in the computation. Other specifications are the same as 
those used in obtaining Fig, 5*2. In this case, however, 
greater frequency intervals at every 50 Hz were used to save 
the computer time instead of 10 Hz interval in Fig, 5*2, 
Although the details are missing in this figure as compared 
to Fig, 5*2, the major effect of the random tension, which 
raises the spectrum values at some frequencies by several 
orders of magnitude, can be seen clearly. The comparison 
with Fig, 5*2 also shows that the frozen-pattern assumption 


200 


frequency, Hz 

Pig. 5.2 Structural displacement spectrum imder the 
frozen-pattern turbulence excitation 



in /Hz 
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is Tmconservative, resulting in gross underestimates in the 
structural response calculation, 

5 A Conclusion 

The turbulence decomposition scheme was applied to the 
two-dimensional problem of a membrane with random tension 
field. The Levy’s series representation was used to describe 
the variation of the wave -number response function, H, in the 
y-direction. Therefore, integrations in the y-direction 
needed to calculate the cross-spectrum of the structural 
response could be carried out simply with only a single inte- 
gration remaining to be done numerically on a computer. The 
effect of the randomness in the structural properties was 
also Investigated, in particular the random non-uniform tension 
in the membrane. It has been shown tlxat with a coefficient 
of variation of only 10^ in the tension field the structural 
response spectrum may be increased by several orders of 
magnitude at some frequencies. The variation of the tension 
field need not be accidental, but it may be caused intention- 
ally to create a spectral pealc, for example, beneficial to 
sicin-friction reduction. This proposition may be investi- 
gated in the future. The suitability of the frozen-pattern 
assumption was, again, examined in this chapter and the same 
conclusion as in Chapter III was drawn; namely, it leads to 
underestimation of structtiral response and omission of some 
important peaks in the spectral density of the response. Thus, 
the use of the frozen-pattern assumption should be avoided for 



72 


the analysis of the structural response spectra under the 
boundary- layer turbulence excitation. 

In Chapter III-Y, the emphasis has been placed on 
structural motions. In the next chapter our focus will be 
shifted to the slcin- friction drag reduction itself. 
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VI, REIHOLDS STRESS ON A COMPLIANT SURFACE 
6,1 Introduction 

From observing the swimming of dolphins and amazed "by 
the smooth motion of the fish Kramer hypothesized that 
favorable interactions between the flabby skin of the fish 
and water could reduce the skin-friction drag on the fish. 

Since then many researches, experimental [36-371! as well as 
theoretical [38-41], have been carried out in this field to 
find the mechanism of this highly applicable phenomenon. 
Aircraft and ship designers are especially interested in this 
problem since reduction of skin- friction will result in a 
decrease in fuel consumption. It has been estimated that the 
skin-friction drag of an aircraft can be as high as 50 ^ of "the 
total drag. Although some experimental studies performed on 
very flexible membranes have indicated friction-drag reduction 
in fully turbulent boundary layers, the theoretical follow-up 
has not been as successful. For historical reviews of previous 
works the reader is referred to [ 36 , 37 ]« 

One requirement for a successful theoretical analysis 
is a thorough understanding of the complicated structural 
motion under the excitation of the boundary- layer turbulence 
[37] » to which Chapters III-V of this thesis were directed. 
Another requirement is the knowledge of the changing flow 
field resulting from the structural motion. In this chapter 
a perturbation approach similar to that used by Ffowcs 
Williams [48] and Slick [4l] will be applied to compute the 
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perturbation Reynolds stress near the structure -fluid interface. 
The analysis will he related directly to the two structural 
models considered in previous chapters. The first model, a 
one-dimensional unsupported beam will be used again to present 
the basic concepts, and then these concepts will be applied to 
the more realistic model of a one-dimensional infinite beam on 
evenly spaced supports. In both cases the turbulence decompo- 
sition scheme will be utilised, 

6,2 Theory 


Pig, 6,1 shows a one -dimensional infinite beam which 
is basically the same beam considered in Chapter III, As the 
beam responds to excitations its motion will generate addi- 
tional pressures in the fluid media on the upper and lower 
sides. As before, denoting these induced pressures by and 
Pgi respectively, the governing equation of the beam motion 
is given by 


A L 

3w 3 w 

m + 11 -^ + El — ~ 

au2 ' at ax'f 


- T 


ax^ 


~ p + (p. - 


P2^z=0 


( 6 , 1 ) 


Equation (6,1) differs from Eq. (3.1) in that a pre-tension T 
and a viscous damping coefficient tj have been included. 

The turbulence pressure field p can be expressed in 
the form previously introduced in Chapter II: 

^i(ukt - kx)' 

A 

— OO 


p(x,t) = 


dP(fc,u) dC(u) 


( 6 , 2 ) 


p (turbulent pressure) 



' ///// //V /////// ///// /rr^T 


Fig. 6.1 One-dijaensioml infinite beam under tbe 
excitation of boundary-layer turbulence 
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She meaning of the functions P and G has been discussed in 
Chapter II, We recall that in the integrand the circular 
freq^uency o) is related to the wave-numher k and the component 
convection velocity u by the relationship to = }cu. Now, the 


induced quantities w, 


Px 


and Pg can be expressed as follows* 


w(x,t) = H(k,u) e 


i(ukt - 


dP(k,u) dG(u) (6.3) 


oo 

P^Cx,y,t) = I P^(k,u,y) " ^^dF(k,u) dG(u) 

(6,/f) 

oo 

PgCxty.t) = I giCulrt - dG(u) 

( 6 . 5 ) 


v/here H, P^, and P^ may be called the wave-number response 
functions for the structural displacement, the radiated 
pressure, and the transmitted sound, respectively. 

To determine the radiated pressure on the upper side 
of the beam, p^, we use the usual assumption that the addi- 
tional pressure can be obtained without regard to the presence 
of the turbulent pressure p. Then, p^^ can be obtained in the 
same way as in Chapter III; 

P^(k,u,0) == 

where TJ^o denotes the free-stream velocity of the flow on the 
upper side of the beam. For the case of incompressible fluid, 
we let 



Then Eq, (6,6) reduces to 


11 


P^(k,.u,0) = p^.k(u - H(k,u) (6.7) 

)?ollowing the same procedure as in Chapter III| the pressure 
generated on the lower side of the "beam is found to he 

« cot(/u^ - ao^ kd/a2) 

F^(k,u,0) = rz 2 H(k,u) (6,8) 

y XI " ^2 ^'^^2 

where the hottom wall of the cavity is assumed to he acousti- 
cally hard, Taking the limit a^^oo, we have the incompress- 
ible flow solution 

Pg(k,u,o) = ° — H(k,u) (6.9) 

k 

For what follov/s v/e shall restrict our attention to the in- 
compressible case. 

Substituting Eqs;* {6,2)-(6,5) » 7) » and (6,9) into 
Eq, (6,1), one obtains 

H(k,u) = £-mu? 4- iTjco + Dk^ t Tk^ - pih.(u - Uoq) 

This equation is similar to those derived in 

however, in these references the effects of the radiated 

pressure p^^ was not taken into account. 

Let the velocity components of the turbulent boundary- 
layer flow on the upper side of the compliant plate be It t u* 
and V + V*, where "u and T are the velocity components when 
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•the boxmdary is rigid and u* and v* are perturbation velocities 
induced by the structural motion. Now assume that the pertur- 
bation velocities, u* and v*, can be expressed as power iSeries 
in yj 

u»(x,y,t) - 2 a (x,t) y^ '(6,11) 

n=0 ^ 

v*(x,y,t) - Z b (x,t) y^ (6,12) 

n=0 ^ 


From the conservation of mass of an incompressible fluid, 
we have 


au* Sv* 
ax','*' ajT 


( 6 . 13 ) 


Tn order that this condition is satisfied the coefficients 
in Eq.s, (6.11) and (6.12) must be related as follows* 


1 

"tJ-Cxit) for n = 1,2, • • * (6,l4) 

n n ax 


Therefore, substituting Eq., {6, Ik) into Eq, (6.12), we have 


v*(x,y,t) = bfv(x,t) - S 


~ 1 3an_i(x,1:) j, 


n=l ^ ax 


(6.15) 


The boundary conditions at the beam surface are 


u(w) t u*(x,w,t) = 0 


( 6 . 16 ) 


v(w) + V* (x,w,t) = 


aw(x,t) 

at 


( 6 . 17 ) 


Taking Tsiylor expansions of u(w) and v(w) about y = 0, and 


79 


neglecting the high order terms of w, one obtains 


u(w) = u{0) + (|”)y,oW + • • . « c||-)y=oW (6.18) 

v(w) =: v(0) + (-|“)y_jjW ■£“••• f» 0 (6,19) 

since u(0) and ^(0) are zero and the continuity of the 
unperturbed fluid requires that 

^3y^y=0 " ^3x^y=0 " ^ 

Substituting Sqs. (6,11) and (6,18) into Eq. (6,l6)j one 
obtains 

2 3u 

a (x,t) + a (x,t)w + a (x,t)w + ' * * = »{-^) w 

ox d o y yziQ 

( 6 , 20 ) 


Keeping only the first terra on the left-hand side of this 
equation, we have 

a^(x,t) = -Uw (6,21) 

where U is defined as 




( 6 . 22 ) 


In a similar way we obtain from Eqs, (6,15)» (6,17) (6,19), 

9w 


aa^CXft) I 3ai(x,t) 
b^(x,t) w - — — r - 


« 1 • — 


3x 


3x 


3t 

(6.23) 


Again, retaining only the first term on the left-hand side, 


bo(x,t) 


3w 

at 


i6.2k) 


Eo[Tiation (6.2^) confirms what might have "been concluded from 
intuitive reasoning that the perturbed velocity in the y- 
dirsction is eq.ual to the structural motion in that direction. 
Substituting Eq.s, (6,21) and (6,24) into Eqs (6,11) and (6,15) » 
we have 


u*(x,y,t) 

v*(x,y,t) 


= - Uvr(x,t) t S a (x,t)y^ 
n=l ^ 

3w(x,t) 3w(x,t) 


at 


+ u 


ax 


n-2 n ax 


n 

y 


(6.25) 


( 6 , 26 ) 


in addition to the continuity eq.uation and the boundary 
conditions, the velocity components in the flow field must 
satisfy two momentum equations, one in the x-direction and 
the other in the y-direction. However, it is sometimes more 
convenient to replace one of these two equations by a diffusion 
eqiiation governing* the vortiqity perturbation, 

av* 

co’(x,y,t) - ^ (6.27) 

In the proximity of the structure this diffusion equation 
is given by 


2 


2 
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where v represents the coefficient of viscosity. We shall 
apply this diffusion equation next to calculate the higher 
order coefficients in the power series of u* and Since 
(O' results from the structural motion it can he expressed 
in a form similar to Bqs, (6,2)-C6,5)» 


(jo*(x,y,t) 



li(k,u,y) 


i(u3ct - }cs) 
e 


dp(k,u) dG(u) 


(6.29) 


Substituting this expression into Eq, (6,28), we have 

2 

— - (k + -^)0(k,u,y) = 0 ( 6 . 30 ) 

ay 

with the condition that the vorticity vanishes at infinity, 

A solution for this differential equation may he expressed as 


_ 2 h 

n(k,u,y) =Q(k,u) exp[-(k + y] ( 6 , 31 ) 


where the function Cl is to he determined. Expanding the 
exponential function on the right-hand side of Eq. (6,31) in 
a 'Jaylor series about y=0, one obtains 


2 

expL"(k + -;y“) 


y] = ? ^ (-1)'' 




iuli 


Ik 

.2 


Therefore, Eq, (6,29) can be written as 




oo 


n 


uj'Cx.y.t) = riridc.u) s -Jr (-d (k‘+^) 


n 


n 


iuk ^2 


JJ 


n=o 


n! 


i(ukt - kx) 


dF(k,u) dG(u) 


(6.32) 


However, from Eqs . (6,25)-(6»27) 


o)*(x»y*t) 


at ax 


a (x,t) 
1 


3x2 


- 2a2(x,t)]y + • • • 


(6.33) 


In consistence with the representations of w and to’, Eq.s« 
(6,3) and (o»29), each also can he expressed as follows* 


Oq 

= J] 


i(ulct - hx) 


dFUc.u) dG(u) (6.3^) 


Substituting this equation and Eq, (6*3) into Eq. (6,33) y 


ua'(x,y,t) = [ul£ H(k,u) - A. (k,u) - Dk H(k,u)y 


- 2A2(k,u)yI! dF(k,u) dG(u) 

. + higher order terms in y (6,35) 

By equating coefficients of y^ in Eqs. (6,32) an^ (6,35), 
we obtain 


n(k,u) 

= uk^H(k,u) - A^(k,u) 

(6,36) 

ii(k,u) 

(k^ + = Uk^H(k,u) 4- 2A2(fe,u) 

(6,37) 


Eliminating XI from these two equations, we have 


J 


I 


S3 


+ 2A2<k,u)] 


(6,38) 


To calculate and Ag we use tiie momentum equation 
in the x-direction 


a(u + B') ,a{u + ii') a(u + u') 

- + (u + u') + (v H* V*) — 


at 


ax 




1 S(p + P, ) J2,_ 

i— + VV (U + U') 

Pi Sx 


(6,39) 


However, for the rigid wall case, we have 


au _ au . _ au i ap ^ 

— — + u — t V - — = - ~ — + "OV u 

at ax By ax 




Thus, by subtracting Eq. (6,^0) from Eq, (6,39) aud neglecting 
the square terms of perturbation velocities, one obtains 


Bu* 

. — Bu* 


, — Bu^ 

. au 


+ u ■“ 

+ u' — 

+ V 


at 

Bx 

ax 

ay 

ay 



Pi 


2 

w u’ 


(6. in) 


It has been shown in Eqs^ (6.18) and (6,19) that at the 
proximity of the beam u and v can be approximated as follov/st 


u(y.) -u(0)> 

_ 

v(y) = v{0) + (^)y=o y + ' • ’ = 0 


Thus, the Wavier-Stokes equation for the perturbation 
velocities near the compliant svrface can be approximated as 



Bli- 




3u* 3u* 1 2 

at" + aF ° 3F 

Substituting Eqs, (6.4), (6.2^), (6.26), and (6.34) in Eq.. 
(6, ,42) and letting y *== 0, we obtain 


ilcP^(Ic,u, 0 ) = -^[uA(k,u) + 2 A 2 (k,u)] ( 6 . 43 ) 

Therefore 

^ Tr TT 2 

A 2 (h,u) = - “ 2 ^ H(}c,u) ( 6 . 44 ) 

Eliminating Ao(h|U) from Eq.s, ( 6 . 38 ) and ( 6 . 44 ), 

At 

A^(lc,u) = uk^H 0 c,u) + ^^)“^P^(lc,u, 0 ) ( 6 . 45 ) 


Now, substituting Eqs. (6.3)» (6.44), and (6,45) into Eos. 
(6,25) and (6,26), we have 


00 

f 


u*(x,y,t) - 


i(ukt - kx) 

B(k,u,y) s dp(k,u) dG(u) 


j 
—00 


(6.46) 


v(x,y,t) = 


i(ukt - kx) 

C(k,u,y) e dP(k,u) dG(u) 


(6,47) 


where 


B(k,u,y) = -UH(k,u) + £uk^H{k,u) + ^ 


fl 



iuk 

+ ) ^ 

V ^ 


P^(k,u,0)}y 4- higher order terms in y (6,48) 
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C(lc,Ujy) = iukH(]£,u) - ikUH(lc,u)y 

+ higher order terms in y (6,49) 

The Reynolds stress associated with the perturbation 
velocities, u* and v', is defined as^ 

rCxiy) = Pj^eCu'v*} (6,50) 

The usual minus sign associated with the definition of the 
Reynolds stress is dropped here since, as shown in Fig, 6,1, 
the positive v* is taken to he opposite to the conventional 
direction. Now, substituting Eq,s. (6,46) through (6,49) into 
Sq., (6,50) and applying Eq. (2.26), one obtains 



# In theory, the imaginary parts in the expressions for u» 
and V* should be aero; however, some approximations have 
been used in the analysis which my lead to non-zero 
imaginary parts in these expressions,^ Such superfluous 
terms should be discarded in calculating the Reynolds 
stress, Therefore, instead of Eq, (6,50) 0^® write 

T(x,y) = Pj^ELHs(u*) Re(v')] 
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pressure, Sp, can "be expressed in terras of the ft^equency 
spectrum using Eq., (2,33)t (6 #51) can he changed also to 

oo 

•cCy) - j| o?(k,Y*y) lp(0,w) W(k - dk dco (6,53) 

— oo 

where 5p(0,o}) denotes the frequency spectrum of the turbulence 
pressure and W accounts for the decay property with respect to 
spatial separation. Due to statistica.1 statiomlity and homo- 
geneity of the structural motion, the first three terms in the 
right-hand side of Eq, (6,53) vanish. Thus the expression for 
the Reynolds stress, reduces to 


T(y) = 


¥p(0,u5) ¥(k - f-) dt dw 


{6.5k) 


We now define a Reynolds number, R, based on the 
wave convection velocity oj/lc and the wavelength 2Tr/kj 


^ ^ (ujA) ( 27T/k) ^ _w_ 

27TV " ^^2 


(6.55) 


Assuming this Reynolds number, R, is large, the square-root 
in the integrand of Eq. {6.5k) is expanded in a Taylor series 
about J /r = Oj 



itth-i _ 









• • ) 


{6.56) 


For values of R much greater than unity, the first term of the 
series alone is adequate. Thus, 
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oo 
— Oo 

^„(0,1U) fOc - f-) dJc du3 (6.57) 

^ Q 

This expression of the Reynolds stress differs from Eq., (33) 
in C^8j in that the radiation pressure, P-|^, takes place of the 
pressure fluctuation. As shomi in the process of obtaining 
^Jq,« (6.^2), the radiation pressure is more directly related 
to the pertin'bation velocities and# therefore, to the perttir- 
bation Reynolds stress than the unperturbed turbulence pres- 
sure on the rigid wall. Thus, this expression of the Reynolds 
stress seems to be more plausible than that in [^^1,^83, 

Now, substitutions Eqs, (6,7) and (6,10) into Eq, 
(6,57), and keeping only the real part which represents the 
physical Reynolds stress, we have 


c 


Tk 


-oo 


- - u.,) - 


2 coth(kd) 2 22 

— ] + ^ } y 


©To(Osfo) dk duo (6,5^) 

It is clear from Eq, (6,5^) that the perturbation Reynolds 

stress vanishes when there is no pressure radiation; i,e«, 

P-(k,— ,0) - G, This happens when the wave propagation speed 
k 

of the structural motion, to/fc, becomes equal to the free-stream 
velocity of the flow, Of course, in the absence of struc- 

tural motion the perturbation Reynolds stress also must be zero. 
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If the beam motion is favorable the Reynolds stress may become 
negative. It has been suggested C^8]| that such a negative 
Reynolds stress in the vicinity of the surface may deprive 
the turbulence in the boundary layer with energy supply and, 
therefore, the turbulence level may decrease. As seen in Eq., 
(6.58) the sign of the perturbation Reynolds stress can be 
altered by changing the structural properties. However, the 
relationship is subtle and extensive numerical studies would 
be required to reach any quantitative conclusions. 


6,3 Periodically Supported Beam 


^Dhe structural model used in the preceding section 
for a preliminary investigation of the change in the Reynolds 
stress in the fluid due to structural motion was an unsupported 
infinite beam. However the results obtained there can be 
extended easily to the case of periodically supported infinite 
beam. This model is shown in Fig, 6.2, and it resembles more 
realistically the construction of an airplane fuselage. 

The governing equation of the beam motion not directly 
on the support is given by 


2 

m “2 “*■ TT “*■ ^ 


O W 


at 


ax 


5 - 


^2 
a w 

ax^ 


= p + (p 

1 



3F=0 


(6,59) 


p (-fcurbuXent prebstira) 



u, u* 


V, V' 


p- (radiated pressure) 
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Pig. 6.2 One-dimensional periodic beam under the 
excitation of boundary-layer turbulence 
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where denotes the viscous damping coefficient! g is the 

loss factor of the beam material, D is the beam rigidity, and 

m, T, p and p are the same symbols used in the preceding 
X ^ 

section. For the purpose of determining the v/ave**number 
response function, H(x,lc,u), the turbulent pressure p should 
be replaced by expCi(ukt - lex)] and the structural response 
w equated with H(x,lc,u) exp(iulct). Other induced quantities 
should also be replaced by their wave-number response functions 
in the same manner. 

Since the supports give rise to multiple reflections 
of the propagation wave in the structure, the same expression 
for the wave-number response function H as that used in Chapter 
IV is suitable; i.e,, 

H(x,k,u) = 21 ■ H exp(-iu x/fi) (6,60) 


where 

+ 2nn = hi + 2 n 7 T (6,61) 

The radiated pressure is governed by the equation 

^2 ^2 

(ft w Pi - Pi = ° 

and subject to the conditions that p^ can propagate only in 
the region y > 0, and ^ that 


^P-t 
/ — ±\ 
^3y 


PlClku ^ It, 1^)2 H 


( 6 , 63 ) 



!Che solution of this system, when evaluated at y=0, is given 


"n 


(tijl - u„)' 


i{uht - 


Tn 


TT n 


(6,64) 


where 


"n = ’^4'n 


(6,65) 


The solution for the case of incompressihle flow can be 
obtained by tahing the limit a^j^-3»o^j i,e. 




p. = S ' e 

i n=-oo 


i(ufet - y:^/X) 


(6v66) 


where 


Pin = Pi#) “n 


( 6 , 6 ?) 


The transmitted pressure Pg in the fluid medium 
•d s y •< 0 is governed by " 


•jS 2 £ 

at^ a^r ^ 


( 6 , 68 ) 


and subject to the conditions 


5y 


= 0 at y “ -d 


aj- = -P2('*) H e 


at y - 0 


( 6 , 6 §) 


( 6 , 70 ) 


The solution for Pg# when evaluated at y = 0, is given hy 


Pg = PgCuk)^ J 


n=-w 


Hn 


(Ukt - /4„s/i ) (g . ) 

Tn 


where 

= (|£)2 . (&)2 ( 6 . 72 ) 

Talcing the limit a^-^ oo, Eq* reduces to the incom- 

pressihle flow solution 


Pg = 


2 oo 
2 

n,=«oo 


oo-thQ^nd/t) i(u!rt - }^„x/I) 


(6,73) 


At this point it is of interest to note that the 
are not all independent of each other due to the constraints 
at the periodic simple supports. At these supports the 
deflections are zero. Thus hy substituting x - 0 into Eq, 
( 6 . 60 ), 

H(0,lc,u) = S = 0 

Xl^-oa ^ 


or 



oo 

- s 


n=t-oo 


n#0 



(6.7^) 


Hence, an alternative representation for H is 


H(x,lc,u) 


n=-oo “ 

n?^0 


e-V*oxA ) 


(6,73) 
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which is also a hetter representation from a computational 
point of view* Substituting Sq.s, (6>66), (6.73) and (6.73) 
into Eq. (6.59) » one obtains 


S Ccp(n) 
n--** 
n#0 


- <p(0) e 






iulct 

e 


_ i(ulct - }^o^/Jt) 


(6.76) 


where 


cp(n) = -m(u!c)^ + i')j(uh) + D(1 + ig)(^)^ + 

^ )VJl‘ 


(6.77) 


To determine the coefficients H we apply the well-lmown 
virtual work principle. Specifically, we assume a virtual 
displacement 


Sw = SH, : e-""^ 

d 


(6.78) 


The virtual work done by the internal and external forces 
must sum up to zero. Bue to the spatial periodicity of the 
structural response and the virtual displacement, it is only 
necessary to apply the virtual work principle to one periodic 
unit. Thus, multiplying Sw to both sides of Eq. (6.76) and 
integrating over 0 ■«£ x c j 2, we have 


cp(3) + cp(0) 2 = -1 for 5 0 

•) IXswao 

n#0 


i^79) 


9^. 


Eq.uation (6,79) shows that the product cp(j) is independent 

tl 

of U i,e,. 


cp(i)J^ j= Cp(2)H2 - * • • = cp(j)Hj = « •. : 

Thus, substituting 

Kjj - K^cp(3*)/cp(n) (6-,8G) 

into Eq, (6,79)* we obtain an equation involving only one 
unlmown 

H^Cp(5)cp(0)^S co^n) “ 3^0 (6*81) 


which is solved readly to give 

Hj = -C(p(3)(p(0)J for a # 0 (6.82) 


Denote the velocity components of the turbulent 
boundary-layer flow on the upper side of the beam again by 
u + u* and v + v" where u and v are the components correspond- 
ing to a rigid wall. In view of Eq, (6,60), the perturbation 
velocities, u* and v*, can be expressed as power series in y 
as follows c 


u"(x,y,t) = S 5 A .y^ e 
n-0 3--00 nj 


i(uht - 


QO OO 


v*(x,y,t) = S .S B .y^ 
n=0 


( 6 . 83 ) 

i6M) 


From Eqs. (6.83) and (6>84-) and the continuity eg.uation for 


9$. 


an incompressible fluids 


Bu* By* 
Bx By 


0 


(4^85) 


one obtains 


„io ^3 ® 


i(t*t - fi.x/i) 


“ n-1 

= S S ny H 
n=l j=-oo 


hd 


i(ukt - ii^xA) n-1 

j ^ ny 


Multiplying both sides of the equation by expd^^x/?.) and 
integrating over 0 < x •< JZ* we have 


BnS Vl,3 ' for n = 1,2, 


(6.86) 


Therefors, by subs'ti'fcu'ting Etj, (6,86) iti*fco Eq, (6,8^), 


w 1 i Ti*-i 

v'(x,y,t) = CBoj + n (-j^) 3 

i(uht “ 


(6-87) 


51he non-slip condition on the plate surface can be 
stated as 


u(w) + u*Cx»w,t) = 0 


C6.8S) 


and 


^ . . Bw(xjt) 

V(w) tV*(X|W,t) = — ^ 


A6.89) 
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Taking Taylor expansions of u(w) and v(w) about y = 0, and 
neglecting the high order terms of w, we have 


(6,91) 

Substituting Eqts. (6.83) and ( 6 i^- 90 ) into Eq.. (6.'88), one . 
obtains 


H(«) = W 

v(w) = 0 


j= 



+ 


+ AajVi* 




(il.) 

'3y^y=0 


OO 



^i(ukt - 


Keeping only the first term on the left-hand side of this 
eq^i^ationi we have 


where 

^ “ «ly5y=0 

Similarlya from Eq.s, (6.87) f (6,89) and (6,91),. 


’( 6 . 92 ) 


(6.93) 


= iukH^ (6.94') 

Now, substituting Eq.s. ( 6 , 92 and (6,94) into Eqis. ( 6 . 83 ) 
and ( 6 , 87 ) » we obtain 


97; 




.S + S A , 

3=-" ' ** n-1 “3 


( 6 . 95 ) 


eo ^M-.! ^ 1 li 

v'Cx.yst) = .E (iulcHi - UH^y -J- S -r -4-3. A .* .y ) 
3 =r_j 8 o 3^3*" n=2 n n-1,3 


iiukt - /t^x/1) 


(6,96) 


The vortici'fcy perturbation 


av* 3u' 

w'(x.y,t) = aJ- - 57 


(. 6 . 97 ) 


must satisfy the diffusion equation near the beams 


Sa>* ,a ^ . 

“ t(— 5 + — g)co' 

at ax"^ ay^ 


(6.58) 


In view of Eqs, (6,55) and (6,96), we may assume that the 
vorticity perturbation has the form 


to'(x,y,t) = .2 n. e 
3=.oo 3 


i(u3ct - p^x/1) 


(6.99) 


Substituting Sq., (6,99) into Eq. '(6,98). multiplying both 
sides of the equation by exp(ijJijx/Ji) , and integrating over 
0 < X ■<J1, one obtains 


iifl. . i:(^)" +^]Q. ^ 0 


(6:^100) 


Under the condition that the vorticity can diffuse only in 
the positive y«domain, the solution of Eq, (6,i00) is found 


to te 


exp£-C(^)^ + ^]“yl 

Expanding the exponential function in aj laylor series about 
y = 0, and substituing the result into Eq* C.6', 99) » we 'have 

^i{ukt - /Jjx/X) 


u>*<x*y,t) ^ S XI 

^TS»0O ' 




However, from the definition of the vorticity, lo* also can be 
•expressed as follows, using Eqs, (J',25) a3ad ( 6 . 96 )^ - 

co’(x,y,t) =: ;S . 

+ higher order terms in y}; • }x^x/SL) 

(6402) 

Equating coefficients of in Eqs* (6.101) and (6.102), one 
obtains 

= Ulc(^)Hj - Ajj (6.103) 

■*■ 2^23 (64ofr) 


Eliminating il. from these equations, we have 

Mj! P^a 2 iulc P '1 ^ 

A^^ = tUc(^)Hj- - C(^) + UHj + 2A2j3 


(6.103) 


Jtoier-Sliokes 

veiociiJies in- tlie vicinity/jof 1;iie *beam can lie bitten aa ; 


SU * ,T ^ 


+ V 


^U' 


(6*100 


Sgs* (6 "60 ^ (6 *^5)^-^ 
and letting y ^ 


lOo- 


. + PjV 




- 0 


(6.io7)- 


Multiplying both sides of the eq.uation by exp(iji^x/C.) and 
integrating over 0 < x <^ l*_one- obtains- 

^23 = - W t’ ^13 - f ^>^^3 

Substitution of Eg.. (6'«1C8) into Eg, (6*105) yields 

A, . = U1C(^)KJ + j 

NoWs substituting Egs. (6,108) and (6,109) into Egs, (6,95) . 
and ( 6 , 96 ), we have 


« ^ i(u3£t - 


( 2 C,y,t) ~ 2 Kj; e 
,i=-oo d 


v»(x,y,t) = .2 S. e’ 
0“-« ^ 


(6;ilG) 


(6,111) 


where 


R3 = -mj + 


■^3 


,2 . 


A.- 


OK,. + J 


3 fi-f'-'X 


I *^3 


iF 


+ higher order terms in y 


i6,itzy 


iu, 

S - = iukHi - TJH^y + hi^er order terms in y (6*113) 
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We note that the velocity components » u* and v*, obtained 
here are induced by a frozen»pattern component j expCi(uht - 

, CO 

hx)3j of the turbulent pressure pj therefore j .S R. exp{ 

oo 3 =-" 3 

and ^ expC-^jz/t ) are the wave-number response 
functions of the velocity components, u* and v% respectively > 
Thus, the total velocities induced by the turbulent pressure 
p-are obtained by superposition as follows* 
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The perturbation Reynolds stress near the beam can 
now be computed from 

T(x,y) “ f^ECu*v*3 (6.116) 

Substituting Eas, (6.1i2)-(6.1i5) into Eg,. (6.116), one obtains 
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As in the pE’eceding section, the integration variahle u in 

( 6 . 117 ) may he replaced hy a)/fc and the wave-nimher spectrnm 
of the turhulence presstire^ Sp, replaced hy the fre<iiiency 
spectrm to glTre 
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vfhere the specific esepression for obtained hy substituting 

Bq.s, (6.112) and ( 6 . 113 ) imto Sq., (6.118) is given as followsi 
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+ higher order terms in y * (6.120) 

Since the structural motion is a stationary random process in 
time, the first term in the expression of. does not con- 
tribute to the integral when substituted into Eq.. (6oll9)* In 
this periodic beam ease, however, the spatial homogeneity of 


the steucttrcal mo*fcioxi Is destroyed 'bseaviae of the supports # 
Consegiiently, the perturhation Reynolds stress near a periodic 
heam is a function of the spatial variable x, A reasonable 
measure of the net effect of the perttirbation Reynolds stress 
appears to be the spatial average defined as follows* 
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Againj define a Reynolds number 


R = 




(6ai22) 


E<iS» C6ell9) and (6*120) may be substituted into Eq., (6*121),,. 
and the result is simplified by e 3 cpanding the square -root 
tetm. in a faylor series about 1 /r = 0, Keeping only the real 
part of the expression which represents the physical average 
Reynolds stress, one obtains 
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Finally, substitution of Eq,s. (6.6?) into Eq,* (6,122) yields 


<r(y)> - 


to 




% /n'g.coJ! „,2j ,2 

n2.« ^2^ -I % ■ “ ' 


5 (O.tu) ¥(lc - dk dco 

P 


(6,12ii.) 


103 


6.^ Concludiijg Hemar&s 

Espressions for the perturbation Reynolds stress 
induced by fluid-solid interaction were obtained for two 
structural models, The first structural model, an infinite 
unsupported beam, was used to develops the basic concepts of 
the analysis. The results were then extended to the second 
model of an infinite beam, simply-supported at equal intervals. 
The second model is a more realistic model for the typical 
fuselage construction of an airplane. The effect of the radi- 
ation pressure was included in the formulation. The pertur- 
bation Reynolds stress involves a double integration with an 
integrand depending on the beam motion and radiated pressure. 
The complexity of the expression permits only qualitative 
discussions in this thesia. Hopefully, quantitative results 
can be documented after extensive numerical studies in the 
future. 

The radiation pressure, p^^, required in the present 
analysis has been obtained from a \mve equation for inviscid 
fluid. More rigorously the viscosity in the fluid and the 
velocity profile in the boundary layer shotild be taken into 
account when determining the radiation pressure. However, at 
the present time no closed form solution is known of this more 
accurate wave equation. 


lo4‘ 

VII. GEMRAL GGNCLUSIOK 

This investigation has "been concerned with the inter- 
action between a turbulent flow and certain types of structinre 
responding to its excitation. The turbulence is typical of 
those associated with a boundary layer, liaving a cross-spectral 
density indicative of convection and statistical decay. It ‘ 
has been shown that a decaying turbulence can be constructed 
from superposing infinitely many components, each of which is 
convecting as a frozen-pattern at a different velocity. This 
turbulence decomposition scheme reduces greatly the computation 
time by r.educing to one-half the number of integration which 
must be perforaed on a computer, Purthermore, the scheme 
provides a convenient way in which experimentally measured 
cross-spectral density of the turbulent pressure fluctuation 
can be incorporated directly into the computation. 

The results of the structure-turbulence interaction 
were jaresented in terms of the spectral densities of the 
structural response and the perturbation Reynolds stress in 
the fluid at the vicinity of the interface. A number of 
structural models were considered in the investigation. 

Among the one-dimensional models were an unsupported infinite 
beam and a periodically supported infinite beam, . The first 
model was used to develops the basic ideas which were then 
applied to the more realistic second model resembling the 
fuselage construction of an aircraft, Por the two-dimensional 
case -e membrane was used to illustrate the type of 


10 $ 


formulation applica’ble to most two-dimensional structures. 
However B a small random variation in the membrane tension 
was included in the analysis since ideally uniform tension 
never exists in practice. Moreover, the mathematical approach 
used in dealing with random membrane tension can be adapted 
to treat other random structural properties in general. Both 
the one -dimensional and two-dimensional structures mentioned 
above were backed by a space filled with an initially 
quiescent fluid to simulate the acoustic environment when 
the structure forms one side of a cabin of a sea- or air- 
craft. 

It has been found that import=»nt spectral peaks of 
the structural response will not appear if decays in the 
turbulence is neglected in the analysis, Ihus, the usual 
Taylor *s hypothesis of fro.zen-pattern turbulence is uncon- 
servative as far as the assessment of structural reliability 
is concerned, The perturbation Reynolds stress is indicative 
of the change in the skin-friction drag due to structural 
motion. It has been shown that* given the statistical 
information of the boundary-layer turbulent pressure field, 
the perturbation Reynolds stress can be altered by varying 
the structui^al parameters. Therefore* the present study is 
potentially useful for designing flight or marine structures 
to minimise the total skin- friction drag. 
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